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ABSTRACT 


The  transient  testing  technique  for  the  determination  of  heat  trans¬ 
fer  characteristics  of  a  porous  medium  has  been  known  since  1929.  It 
requires  a  good  deal  of  care  and  one  of  the  important  restrictions  is  a 
step  change  in  temperature  must  be  applied  to  the  medium  at  the  beginning 
of  a  test  run.  This  thesis  involves  the  design,  test  and  evaluation  of 
an  experimental  set-up  that  utilizes  an  "integral  heater"  to  accomplish 
the  step  change  in  temperature.  In  addition,  there  is  presented  heat 
transfer  and  flow  friction  characteristics  for  a  20°  skew  matrix  con¬ 
structed  from  perforated  nickel  material.  This  data  is  compared  to  a 
20°  skew  matrix  made  from  stainless  steel  material  that  has  been  pre¬ 
viously  reported  by  C.  P.  Howard.  Also,  data  is  presented  for  two  "paral¬ 
lel  plate"  matrices,  one  constructed  from  brass  sheet,  the  other  from 
perforated  nickel  material.  A  comparison  of  their  relative  performances 
is  made. 
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TNT'IODUCTION. 


1. 

A  compact  heat  exchanger  is  one  not  necessarily  of  small  bulk  and 
weight,  but  one  which  incorporates  a  heat  transfer  surface  having  a  high 
area  density. 

Arbitrarily,  it  has  been  specified  that  ^  —  200  ft^/ft^  makes 

for  a  compact  heat  exchanger,  where  >0  is  the  ratio  of  the  heat  trans- 
2  3  + 

fer  area  to  core  volume,  ft  /ft  .  (15) 

Motivation  for  using  compact  surfaces  is  to  gain  high  exchanger 
effectiveness  in  given  box  volume  and  weight  limitations  and  at  a  reason¬ 
able  cost.  Great  impetus  has  been  given  the  compact  heat  exchanger  field 
by  recent  applications  as  a  regenerator  in  aircraft  turboprop  power  plants, 
as  a  regenerator  in  the  gas  turbine  automobile  engine  and  by  their  use 
in  cryogenic  plant  exchangers  in  the  production  of  liquid  oxygen  and 
nitrogen.  As  the  applications  for  compact  heat  exchangers  increase, 
design  problems  increase.  A  successful  compact  heat  exchanger  design 
and  development  program  requires  access  to,  or  the  ability  to  generate 
knowledge  of  the  basic  heat  transfer  and  friction  behavior  of  compact 
surfaces  of  interest.  It  is  also  desirable  that  the  designer  would  have 
available  as  much  information  on  surface  performances  as  possible  to  re¬ 
main  flexible  in  his  design  with  regard  to  the  limitations  that  have  been 
placed  upon  him. 

Building  a  complete  heat  exchanger  to  evaluate  a  new  surface  as  to 
its  heat  transfer  and  friction  performance  would  be  prohibitively  ex¬ 
pensive.  The  alternate  method  is  to  generate  the  performance  character¬ 
istics  of  a  surface  of  interest  by  testing  small  matrix  modules  and  use 

+Numbers  in  parentheses  refer  to  references  in  the  bibliography. 
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this  data  to  design  a  prototype  heat  exchanger. 

Briefly,  the  objectives  of  this  thesis  are:  to  build  an  experimental 
rig  that  will  test  small  matrix  modules  to  determine  their  heat  transfer 
and  flow  friction  characteristics,  and,  to  test  matrices  made  from  per¬ 
forated  nickel  material  to  determine  these  characteristics. 

One  of  the  simplest  and  most  satisfactory  methods  of  measuring  the 
forced  convection  heat  transfer  characteristics  of  a  matrix  utilizes  a 
transient  test  technique  or  so  called  single  blow  method.  This  test 
technique  consists  of  heating  the  homogeneous  porous  solid  of  knovm 
weight  and  geometry  with  a  uniformly  distributed  heat  source  and  then 
at  a  given  instant,  removing  the  heat  source  and  simultaneously  cooling 
the  solid  with  a  fluid  at  a  uniform  lower  temperature.  The  heating  and 
cooling  steps  can  be  reversed  and  the  same  results  achieved.  A  time- 
temperature  response  curve  of  the  fluid  stream  leaving  the  solid  is 
recorded.  By  comparing  this  record  to  theoretically  developed  non-dimen¬ 
sional  temperature  response  curves,  the  heat  transfer  characteristics  of 
the  solid  can  be  determined. 
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2. 


SUMMARY  OF  THEORY. 


A  theoretical  analysis  of  the  transient  heat  transfer  behavior  of 
a  porous  medium  was  presented  by  Schumann  (2)  in  1929.  A  porous  medium 
may  be  defined  as  an  aggregation  of  solid  bodies  so  arranged  that  the 
void  spaces  between  the  bodies  are  connected  and  will  permit  the  passage 
of  a  fluid.  This  definition  can  of  course  be  extended  to  include  the 
matrix  surfaces  utilized  in  a  compact  heat  exchanger. 

In  Schumann's  analysis,  the  porous  medium  is  considered  to  be  homo¬ 
geneous  and  initially  at  a  uniform  temperature,  t^.  A  fluid  at  the  same 
temperature  is  flowing  through  the  medium.  At  a  certain  instant  the 
temperature  of  the  entering  fluid  is  assumed  to  change  to  a  higher  tem¬ 
perature  t^^.  The  problem  is  then  to  find  the  temperature  of  the  fluid 
and  the  solid  as  functions  of  time  and  position  in  the  matrix  or  porous 
media.  The  basis  of  the  analysis  are  the  theoretical  generalized  heat¬ 
ing  curves  that  are  developed  from  an  energy  balance  on  the  matrix. 


Solid 


3 


Assume ; 


1)  Properties  of  the  fluid  are  independent  of  temperature, 

2)  Fluid  flow  is  steady. 

3)  The  porous  solid  is  homogeneous. 

4)  The  thermal  conductivity  of  both  fluid  and  solid  is  infinite 
perpendicular  to  the  flow  direction. 

5)  Thermal  conductivity  of  fluid  is  zero  in  the  flow  direction. 


Heat  absorbed  by  the  solid  =  ^  4s  41s 

Heat  transferred  to  solid  by  convection  =  h  b  -  "ti)  cl K 

.  I 

Heat  transferred  from  the  fluid  by  convection  =  mCp  "gxr 
Heat  transferred  from  solid  by  conduction  =  "  As 


An  energy  balance  for  the  fluid  and  solid  results  in; 


nCp  ^  Jk  -h  -  ts)  =0 


=  LL 

HCs 

z  =  generalized  position  variable  =  N 


Let;  T=  generalized  time  variable  =  0 


(1) 

(2) 


N 


tu 


h  A 
m  c 


m  c  L 
P 


Making  the  above  substitutions,  equations  1)  and  2)  become; 

:i  (±s-U) 

The  experimental  method  used  requires  that  the  above  assumptions 


and  the  following  boundary  and  initial  conditions  be  closely  observed. 
These  conditions  are  as  follows; 

1)  Initially  the  core  is  at  a  uniform  temperature. 

2)  At  time  t  =  0,  the  temperature  of  the  incoming  fluid  changes 
instantaneously. 
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3)  System  is  adiabatic,  that  is,  only  heat  is  transferred  be¬ 
tween  the  matrix  core  and  the  fluid, none  with  the  surroundings. 

With  the  idealizations,  initial  and  boundary  conditions,  and  the  additional 
assumption  that  the  thermal  conduction  in  the  solid  is  zero  parallel  to 
the  direction  of  fluid  flow  (i.e.,  A  ,  the  material  conduction  para¬ 
meter,  =  0),  Schumann's  solutions  are  as  follows: 


-(z  +T  ) 


CO 

2 


d(zr)" 


(2i  nIt  z) - 


1  - 


+  T) 


oo 


n^o 


n 


z 


d” 

d(zT  )" 


Jo  (2ipr7) 


This  solution  was  first  used  as  the  basis  for  a  transient  technique  by 
Furnas  (16)  in  1932.  The  fluid  temperature  was  measured  at  the  exit  of 
the  matrix  where  x  =  L,  therefore  z  =  The  experimental  data  was 

compared  with  Schumann's  theoretical  constant  z  curves.  The  theoretical 
z  curve  which  best  fit  the  experimental  data  was  considered  as  the 
of  the  matrix  and  from  this  value,  h,  the  convection  conductance  was 
determined.  The  biggest  disadvantage  of  this  method  is  in  the  calcula¬ 
tion  of  an  infinite  series  for  each  point  and  the  number  of  terms  neces¬ 
sary  for  a  good  approximation  increases  as  increases. 

In  1950,  G.  L.  Locke  (1)  developed  another  method  for  the  evaluation 
of  experimental  from  the  theoretical  curves.  In  this  method,  only 

the  maximum  slopes  of  the  theoretical  and  experimental  curves  need  to  be 
compared.  He  derived  an  expression  for  the  slope  of  the  generalized  heat¬ 
ing  curve,  (t^  -  as  a  function  of  Locke's 
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expression  shows  that  the  maximum  slope  of  the  generalized  heating  curve 
is  a  unique  function  of 


d(  f2  ’-j) 

IfLliu 

dTrTSD 


2 

^’tu 


-  J^(2i  )  e"^\u 


^  ) 


The  above  implies  that  the  downstream  fluid  temperature  is  measured  at 

X  =  L.  Therefore,  z  =  N  and  t-  =  t.„.  This  method  has  become  known 
*  tu  f  f2 

as  "The  Maximum  Slope  Method". 

One  of  the  idealizations  made  by  Schumann  and  followed  by  Locke  was 
that  the  thermal  conductivity  of  the  solid  parallel  to  the  flow  is  zero. 
This  is  very  difficult  to  achieve  in  an  experiment  with  a  continuous  flow 
channel  in  the  matrix.  Although  the  effect  of  finite  thermal  conductiv¬ 
ity  parallel  to  the  flow  is  negligible  at  high  Reynold's  numbers,  it  can¬ 
not  be  ignored  at  low  Reynold's  numbers.  Recent  works  by  Crestwick  (3), 
Mondt  (4),  and  Howard  (5)  consider  the  effects  of  longitudinal  thermal 
conduction  in  the  solution  of  Schumann's  transient  heating  problem.  To 
correlate  results,  they  utilize  the  conduction  parameter,  A  .  This 
parameter  is  the  ratio  of  conduction  heat  transfer  to  the  change  in  in¬ 
ternal  thermal  energy  of  the  flowing  fluid  for  the  same  temperature  differ¬ 
ence.  \<Jhen  A  =  0,  the  maximum  slope  is  a  unique  function  of  ^s  put 

forth  by  Locke,  however,  when  X  ^  0,  the  maximum  slope  is  a  function  of 

N  and  A  . 
tu 

Howard,  in  his  paper  (5),  included  graphs  and  a  table  of  maximum 
slopes  as  a  function  of  and  A  .  He  derived  these  values  by  finite 
difference  calculations  on  a  digital  computer  and  proposed  that  the  values 
of  maximum  slopes,  as  presented,  are  within  +  2%  for  a  given  N  and  A  . 
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Ideally  then,  if  the  idealizations,  initial  and  boundary  conditions 
can  be  adhered  to,  all  that  is  required  for  the  experimental  evaluation 


of  the  heat  transfer  characteristics  of  a  matrix  is  to: 

1)  Determine  the  conduction  parameter  A  . 

2)  Experimentally  determine  the  maximum  slope  of  the  heating  or 
cooling  curve. 
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3. 


OBJECTIVES. 


In  the  past,  there  has  been  a  variety  of  expertments  conceived  to 
utilize  the  principles  previously  discussed,  by  approaching  the  idealiza¬ 
tions,  Initial  and  boundary  conditions  as  closely  as  possible.  Perhaps 
the  most  difficult  condition  to  achieve  is;  "at  time  t  =  0  the  tempera¬ 
ture  of  the  incoming  fluid  changes  instantaneously". 

Locke  (1)  and  Coppage  (9)  both  employed  an  air  switching  device  with 
which  air  could  be  quickly  directed  through  the  matrix  or  by-passed  to 
the  atmosphere.  By  use  of  this  switch  they  imposed  a  "step  change"  in 
temperature  on  their  test  matrices.  To  approach  the  "step  change"  in 
temperature,  "requires  a  fast  operating  fluid  switch  for  introducing 
the  hotter  (or  colder)  fluid  in  order  to  assure  an  instantaneous  change 
of  fluid  temperature.  However,  even  with  rapid  switching,  one  of  the 
major  difficulties  here  is  the  Inevitable  thermal  capacity  of  ducting, 
switch  device  and  matrix  holding  device  causes  unavoidable  changes  in  the 
entering  temperature".  (9) 

Solar  Aircraft  Company  recently  did  some  work  testing  core  modules 
for  their  T-600  engine  (10).  In  their  experimental  rig,  they  utilized 
two  steady  streams  of  air  and  a  carriage  tray  that  houses  the  core  module. 
The  tray  was  inserted  through  the  side  of  the  test  cabinet  and  was  moved 
from  a  warm  air  stream  to  a  ambient  or  cold  air  stream  to  simulate  a  step 
flunction  change  of  temperature.  This  type  of  switching  was  done  manu¬ 
ally. 

Howard,  in  his  work  with  new  compact  surfaces  (6),  employed  a  "slid- 
ing  drawer".  This  drawer  holds  the  matrix  and  the  downstream  thermo¬ 
couples.  The  drawer  can  be  inserted  or  withdrawn  from  the  test  section 


through  which  ambient  air  flows  continuously.  The  matrix  is  heated 
when  the  drawer  is  in  the  withdrawn  position  by  an  independent  tempera¬ 
ture  controlled  air  stream. 

One  of  the  objectives  of  this  thesis  was  to  construct  and  test  an 
experimental  rig  that  would  improve  upon  those  previously  used.  This 
test  apparatus  described  herein  differs  from  other  test  facilities 
chiefly  in  the  method  used  to  simulate  the  initial  condition  of  a  step 
change  in  temperature  of  the  fluid  (air). 

The  test  apparatus  constructed  has  but  one  air  stream,  the  air  is 
heated  by  a  high  resistance  nichrome  wire  heater  complex,  an  integral 
part  of  the  apparatus,  located  just  upstream  of  the  test  matrix.  The 
time  constant  or  response  of  the  heaters  is  at  least  an  order  of  magni¬ 
tude  smaller  than  the  response  of  the  thermocouples  and  recording  devices 
used.  By  turning  the  heaters  on  or  off,  the  step  change  in  temperature 
is  simulated.  The  other  systems  required  some  manual  dexterity  on  the 
part  of  the  experimenter  to  begin  a  data  run,  whereas  this  system  in¬ 
volves  only  turning  off  or  on  the  heaters.  The  temperature  recording 
device  is  started  prior  to  the  start  of  a  run.  The  results  achieved  by 
this  system  were  an  improvement  over  other  systems.  With  the  "sliding 
drawer"  and  "quick  changing  valve"  there  is  a  change  in  the  flow  and  there 
fore  the  hydrodynamic  as  well  as  the  thermal  boundary  layer  must  develope. 
In  this  system,  the  hydrodynamic  boundary  layer  is  already  formed,  only 
the  thermal  boundary  layer  is  affected.  A  complete  description  of  equip¬ 
ment  is  presented  in  APPENDIX  B. 

In  addition  to  the  construction  of  this  test  apparatus  and  the  verifi 
cation  of  its  reliability  by  comparing  with  published  data,  extension  of 
the  curves  published  by  Howard  (6)  into  the  high  Reynold's  number  region 
was  attempted. 
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Finally,  heat  transfer  and  flow  friction  characteristics  of  several 


new  surfaces  were  determined. 
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4.  EXrERIMliNTAT.  TECHNIQUES 

In  reference  (5)  Howard  presents  results  from  a  numerical  finite- 
difference  method  of  calculating  the  transient  behavior  of  a  porous 
media  when  subjected  to  a  step  change  in  fluid  temperature.  These  cal¬ 
culations  consider  the  case  where  longitudinal  thermal  heat  conduction 
cannot  be  neglected.  These  data  are  intimately  involved  in  the  determina¬ 
tion  of  and  subsequently  the  Colburn  j  factor  as  presented  in  this 
thesis.  For  this  reason  the  restrictions,  idealizations  and  conditions 
Howard  placed  on  his  solution  will  be  restated  here.  It  is  these  re¬ 
strictions,  idealizations  and  conditions  that  the  experimental  set-up 
must  approach  as  closely  as  possible  to  achieve  reliable  results. 

1)  The  flow  entering  the  matrix  is  steady  and  uniform  In  velocity 
and  temperature,  and  remains  steady  and  uniform  at  any  cross 
section  as  it  flows  through  the  matrix. 

2)  The  thermal  conductivity  of  the  matrix  is  infinite  in  the 
direction  normal  to  the  flow  and  finite  in  the  direction  para¬ 
llel  to  the  flow. 

3)  The  thermal  capacity  of  the  fluid  contained  at  any  time  within 
the  matrix  is  small  compared  with  the  thermal  capacity  of  the 
matrix. 

4)  The  thermal  properties  of  the  fluid  and  matrix  are  constant 
and  uniform. 

5)  The  convective  heat  transfer  coefficient  is  some  suitable 
average  and  remains  uniform  and  constant. 

6)  At  time  zero  the  change  in  fluid  temperature  will  be  a  step 
function,  with  the  matrix  and  its  entrained  fluid  initially  at 
some  uniform  and  equal  temperature. 

The  first  idealization  was  attained  by  careful  experimentation. 
Uniform  velocity  profiles  were  obtained  by  wire  screen  flow  straighteners . 
Uniform  temperature  profiles  were  obtained  by  a  wire  screen  type  heater 
explained  in  detail  in  APPENDIX  B.  The  velocity  profiles  were  taken  at 
the  inlet  of  the  matrix  over  entire  range  of  flow  rates  to  ensure  that 
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the  profiles  were  uniform.  The  temperature  profiles  were  taken  over  the 
entire  range  of  flow  rates  and  found  to  be  uniform  within  +  .5®?. 

The  second  idealization  cannot  be  strictly  realized  in  practice.  The 
effect  of  the  deviations  from  this  idealization  can  be  made  negligible  in 
the  solid  if  the  resistance  to  heat  transfer  within  the  solid  is  small 
compared  to  the  heat  transfer  resistance  between  the  solid  and  fluid. 

The  third  idealization  restricts  the  fluid  to  a  gas  rather  than  a 
liquid. 

The  fourth  idealization  can  be  realized  if  the  temperature  variation 
is  limited  to  20  -  25®F.  rise  above  ambient.  This  would  place  the 
operating  temperature  range  at  60  to  80°F.  In  this  range  the  thermal 
properties  of  the  fluid,  particularly  a  gas,  vary  insignificantly  and  can 
be  considered,  as  a  good  approximation,  constant  and  uniform. 

The  fifth  idealization  is  imposed  upon  the  computer  solution  and,  of 
course,  the  determination  of  the  convective  heat  transfer  coefficient  is 
the  ultimate  result  of  the  experimental  process. 

Finally,  two  sets  of  thermocouples  were  placed  one  upstream  and  one 
downstream  of  the  matrix.  Prior  to  the  beginning  of  each  run,  these 
thermocouple  complexes  were  read.  When  a  steady  temperature  differential 
was  indicated,  the  run  was  begun. 

Herein  lies  the  improvement  over  the  "sliding  drawer"  technique.  A 
series  of  nichrome  wire  heaters  were  installed  in  the  air  stream  to  heat 
the  air  and  subsequently  the  matrix.  To  begin  a  run,  the  heaters  were 
turned  on  or  off,  and  the  temperature  vs.  time  recorder  started.  The  time 
constant  (TC)  of  the  nichrome  wire  heater  is  very  small  (TC  =  .0425  sec 
for  m  =  1000  Ib/hr,  TC  =  .259  sec  for  m  =  10  Ib/hr)  and  hence  the  approach 
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to  the  final  idealization;  at  time  t  =  0  the  change  in  fluid  temperature 
will  be  n  step  function. 

A  schematic  diagram  of  the  experimental  system  in  which  air  is  used 
as  the  working  fluid  is  shown  in  Figure  1.  Figure  2  is  a  photograph  of 
the  entire  test  apparatus.  Essentially,  the  system  consists  of  a  source 
of  air,  an  air  metering  device,  flow  straightener  screen  pack,  an  air 
heater  section,  a  test  matrix,  a  pressure  measuring  system  and  a  tempera¬ 
ture  measuring  system. 

Air  is  drawn  through  the  apparatus  by  a  six-stage  centrifugal  air  com¬ 
pressor.  Air  flow  rates  are  measured  by  an  ASME  standard  D  and  T)/2  ori¬ 
fice  meter  with  changeable  orifice  plates.  Mass  flow  through  the  system 
is  controlled  by  two  sliding  valces,  one  on  the  inlet  to  the  compressor 
and  one  on  the  compressor  outlet.  In  addition,  there  is  a  fine  control 
by-pass  valve  located  just  ahead  of  the  compressor.  As  the  air  passes 
through  the  orifice  plate  it  comes  to  the  screen  flow  straightener  at 
which  point  the  flow  cross-section  changes  from  round  to  square.  The  air 
then  flows  through  the  air  heater  section  which  serves  a  dual  purpose  of 
straightening  the  air  velocity  profile  and  heating  the  air.  In  this 
section  the  air  is  heated  uniformly  to  a  predetermined  temperature,  usually 
20°F  above  ambient.  The  heated  air  passes  into  the  matrix  and  heats  the 
matrix  to  the  same  uniform  temperature.  At  this  point,  a  run  is  ready  to 
start.  By  referencing  the  thermocouples  downstream  of  the  matrix  to  the 
thermocouples  upstream  of  the  heaters,  the  initial  temperature  differen¬ 
tial  can  be  very  closely  controlled.  This  tem.perature  differential  is 
continuously  recorded  by  a  Leeds  and  Northrup  "Speedomax"  strip  chart 
recorder  when  a  run  is  started.  The  recorded  trace  has  a  distinct  ad¬ 
vantage  in  that  there  is  no  transposition  of  data  necessary,  hence,  a 
decrease  in  the  uncertainty  of  the  results.  The  resultant  trace  or  curve 
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is  the  primary  data  for  determining  j.  The  of  the  matrix  can  be  evalu¬ 
ated  by  determining  the  maximum  slope  of  this  fluid  temperature  differ¬ 
ence  curve  during  the  cooling  transient.  This  is  the  basis  of  the  maxi¬ 
mum  slope  technique  as  developed  by  Locke  (1), 

The  flow  friction  performances  of  the  matrices  tested  were  obtained 
by  determining  the  isothermal  pressure  drop  across  the  matrix  and  the 
mass  flow  rate  of  the  fluid  through  the  matrix.  The  pressure  drop  data 
were  obtained  from  static  pressure  taps  located  in  the  test  section  im-  , 
mediately  upstream  and  downstream  of  the  test  section.  The  entering 
fluid  static  pressure  was  also  recorded.  The  mass  flow  rate  data  were 
obtained  in  the  air  metering  section  as  previously  described.  In  summary, 
the  following  data  were  recorded  for  each  run: 

a)  Orifice  diameter  (D  ) 

b)  Pressure  drop  across  the  orifice  (  AP  ) 

c)  Pressure  drop  across  the  matrix  (  AP  ? 

d)  Entering  static  pressure  (Ps)  ™ 

e)  Orifice  temperature  (T  ) 

f)  Atmosphere  pressure  (P  ) 

g)  Temperature-Time  cooling  or  heating  curve 

h)  Temperature  recorder  chart  speed 

A  complete  description  of  the  experimental  apparatus  and  equipment  is 
in  APPENDIX  B.  APPENDIX  A  discusses  the  relation  of  the  data  to  the 
equations  governing  j  &  f. 
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5.  DESCRIPTION  OF  SURFACES. 

There  are  data  presented  in  this  thesis  for  seven  surfaces.  Of 
the  seven,  three  are  new.  Three  of  the  remaining  four  surfaces  were 
tested  by  Howard  (6).  The  remaining  surface  'was  used  primarily  as  a 
test  surface  to  check  the  reliability  of  the  experimental  apparatus. 

The  test  surface  was  a  matrix  consisting  of  20  randomly  stacked 
stainless  steel,  60  x  6C  mesh  screens.  A  photograph  of  the  matrix 
along  with  a  tabulation  of  its  physical  and  geometrical  properties  may 
be  found  in  Figure  3.  The  heat  transfer  and  flow  friction  characteris¬ 
tics  for  such  a  matrix  were  published  by  Coppage  in  1952  9),  and  was 

considered  a  good  check  on  the  performance  of  the  experimental  set-up. 

The  three  matrices  examined  by  Howard  and  re-examined  in  this  thesis 

were : 

1)  Plate-fin  reference  matrix 

2)  10°  skew  passage  matrix 

3)  20°  skew  passage  matrix 

These  particular  surfaces  were  chosen  to  provide  an  additional  check 
on  the  reliability  of  the  experimental  setup  and  to  provide  a  compari¬ 
son  between  the  "integral  heater"  and  "sliding  drawer"  experimental  techni¬ 
ques. 

The  plate-fin  surface  is  formed  by  corrugated  sheets  separated  by  a 
flat  spacing  sheet  to  prevent  meshing  of  the  corrugations.  This  matrix 
was  to  serve  as  a  reference,  since  the  theoretical  laminar  flow  solutions 
for  both  heat  transfer  and  frictional  behavior  of  triangular  cross-sec¬ 
tions  are  available.  However,  difficulty  was  encountered  in  forming  the 
material  into  triangular  cross-sections,  and  therefore  it  fails  as  a  re¬ 
ference,  but  the  results  were  interesting  and  will  be  discussed  later. 
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Figure  4  presents  photographs  of  the  plate-fin  matrix  and  a  tabulation 
of  its'  physical  and  geometrical  properties. 

The  skewed  passage  also  employs  corrugated  sheets,  but  unlike  the 
plate-fin,  the  spacer  plate  is  not  included.  To  prevent  meshing  of  the 
corrugations,  the  sheets  are  skewed  relative  to  one  another  forming  what 
is  defined  as  the  "skew  angle".  A  photograph  of  the  matrix  and  a  tabula¬ 
tion  of  its  properties  may  be  found  in  Figure  5. 

Of  the  three  new  surfaces  tested,  one  is  a  version  of  the  previous¬ 
ly  tested  20°  skew  matrix  using  perforated  or  slotted  nickel  plates  for 
a  material  vice  a  continuous  material.  The  plates  were  corrugated  on 
the  same  device  used  for  the  original  20°  skew  to  maintain  the  same  corru¬ 
gation  geometry.  The  perforations  should  reduce  stagnant  areas  in  the 
matrix,  thereby  increasing  the  heat  transfer  characteristics. 

Finally,  two  matrices  of  a  parallel  plate  design  were  constructed 
and  tested.  The  "parallel  plate"  matrix  is  similar  to  the  plate-fin 
in  that  a  spacer  plate  is  placed  between  the  formed  sheets  to  prevent  mesh¬ 
ing.  The  sheets  are  formed  by  rolling.  When  assembled,  the  matrix  is  a 
maze  of  rectangular  passages  with  an  aspect  ratio  of  77.  Aspect  ratio  is 
defined  in  reference  (8)  as  the  ratio  of  the  two  sides  of  a  rectnagular 
flow  passage,  b/a.  The  sketch  below  is,  ideally,  what  was  intended  in 
designing  the  "parallel  plate"  matrix. 


One  matrix  was  formed  from  a  brass  material  and  is  used  as  a  refer¬ 
ence  to  theoretical  laminar  flow  data  for  parallel  configurations  (8). 

The  second  matrix  was  formed  from  the  perforated  nickel  plate  to  examine 
the  effects  of  the  perforations  on  this  type  of  geometry.  Pictures  and  a 
tabulation  of  properties  for  the  "parallel  plate"  matrices  are  in  Figure  6. 
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6. 


PRESENTATION  OF  RESULTS. 


The  heat  transfer  and  flow  friction  data  for  each  matrix  are  pre¬ 
sented  in  Tables  I  to  VII,  These  data  are  also  presented  in  graphic 

2/3 

form,  the  Colburn  heat  transfer  factor,  j  (  =  )>  friction 

factor,  f,  are  plotted  as  functions  of  Reynold's  Number,  N^,  in  Figures  7 
through  16.  Plotted  also  in  Figures  13  through  15,  for  comparison  pur¬ 
poses,  are  the  theoretical  laminar  flow  solutions  for  parallel  plates. 

For  the  heat  transfer  data,  Prandtl  number  (Np^^)  was  not  really  a 

test  variable,  since  the  working  fluid  was  air  at  moderate  temperatures. 
2/3 

Npp  is  included  as  an  approximation  to  the  effect  of  Prandtl  number 

over  a  moderate  range  of  Prandtl  numbers.  Basic  laminar  boundary  layer 
solutions  indicate  that  in  the  Prandtl  number  range  of  0.5  to  1.0  (gases), 
the  Prandtl  number  enters  the  solution  as  approximately  the  2/3  power. 

Since  a  moderate  temperature  range  is  used  in  the  experiment,  air 
properties  throughout  the  matrix  are  practically  constant.  Therefore, 
all  the  properties  of  air  used  in  the  calculations  of  the  heat  transfer 
and  friction  characteristics  were  those  of  air  at  the  bulk  temperature  of 
the  air  stream  as  measured  at  the  orifice  inlet. 

The  Reynold's  number  is  based  on  the  hydraulic  diameter  and  mass  flow 
velocity. 


A 


The  hydraulic  diameter  is  defined  as: 


D 


H 


=  4r^  =  4( 
=  4 

=  4 


Flow  Passage  Cross-sectional  area  , 

Flow  Passage  Wetted  Perimeter 

Void  Volume  x  _ L _ 

L  Heat  Transfer  Area 
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The  mass  flow  velocity  is  based  on  the  matrix  free  flow  ctoss-sectional 
area : 


G  = 


pA 


The  friction  factor  includes  the  entrance,  exit  and  flow  accelera¬ 
tion  effects.  The  entrance  and  exit  effects  normally  provide  only  a 
small  contribution  to  the  overall  pressure  drop  in  a  matrix.  Core  fric¬ 
tion  is  the  terra  that  controls  the  magnitude  of  the  pressure  drop.  As  a 
result,  high  accuracy  in  the  evaluation  of  the  entrance  and  exit  loss  co¬ 
efficients  is  not  required.  These  coefficients  and  are  functions 
of  flow  geometry  and  Reynold's  number  and  were  obtained  from  Chapter  IV, 
Figures  18  to  21,  as  appropriate,  reference  (8). 
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7.  UNCERTAINTIES  AND  LIMITATIONS. 

As  put  forth  in  the  theory,  the  experimental  method  utilized,  re¬ 
quires  great  care  to  approach  the  idealizations  and  initial  and  boundary 
conditions  of  the  theoretical  analysis.  The  magnitude  of  errors  intro¬ 
duced  if  the  experimental  system  does  not  duplicate  the  ideal  system  is 
difficult  to  evaluate.  Each  of  the  conditions  and  idealizations  were 
carefully  considered  and  the  experimental  system  was  designed  to  approach 
as  closely  as  possible  the  ideal  system.  These  items  are  discussed  in 
the  "Experimental  Technique"  section.  No  attempt  will  be  made  here  to 
evaluate  or  estimate  this  type  of  error. 

Other  experimental  errors  are  caused  by:  1)  uncertainty  in  the 
value  of  physical  constants,  2)  inaccuracy  in  geometric  measurements 
and  3)  instrument  inaccuracies.  The  method  employed  herein  to  describe 
the  uncertainty  in  the  values  reported  is  explained  fully  in  Reference  (14). 
However,  before  stating  these  uncertainties,  an  example  of  the  process 
follows; 

A  theorem  of  statistics  as  stated  in  Reference  (14)  is: 

"If  R  is  a  linear  function  of  n  dependent  variables,  each  of  which 
is  normally  distributed,  then  the  relation  between  the  interval  for  the 
variables,  w^,  and  the  interval  for  the  result,  wj^,  which  gives  the  same 
odds  for  each  of  the  variables  and  for  the  result  is; 


where  w^  =  uncertainty  interval  based  on  certain  odds. 
Consider  the  Colburn  j  factor 


j  =  N, 


^N 

tu  A  PR 


by  the  theorem  and  letting  P  =  N_ 


2/3 
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if  both  sides  are  divided  by  j  then: 


'til  = 

J 


(1) 


Before  an  equation  such  as  this  can  be  evaluated  the  uncertainty 
terms  must  be  known. 

The  physical  constants  necessary  are  listed  below.  References  (12) 
and  (13)  were  sources  for  the  values  of  these  constants.  The  uncertainties 
in  these  values,  as  closely  as  can  be  determined,  are  also  listed  below. 

c  +  0.57o 

s  — 

c  +  0.57o 


+  1.07, 
+  2.07, 


k  +  0.57, 

s  - 


The  uncertainty  in  physical  linear  measurements  could  be  consider¬ 
ed  negligible  if  core  construction  did  not  have  minor  inconsistencies. 

The  uncertainty  in  the  weight  of  a  core  is  considered  negligible.  Conser¬ 
vative  estimates  of  the  uncertainty  in  linear  measurements  are; 

A,  ,  A  ,  A  ,  +  1.07, 

fr’  c’  s’  - 

L  +  0.57. 

W  <  0.17, 
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Instrumentation  errors  were  such  that  only  the  manometers  need  be 
considered.  With  the  exception  of  the  orifice  temperature,  all  tempera¬ 
tures  measured  were  temperature  differences  and  these  were  measured  in 
inches  as  indicated  in  Appendix  A.  Since  pressures  varied  and  manometers 
to  measure  them  changed,  the  maximum  uncertainty  encountered  was  used  in 
the  error  analysis  and  is  listed  below. 

P  +  1.257o 

o  — 

P  +  1.70% 

m  — 

''atr.  negligible 

Another  uncertainty  which  requires  special  mention  is  that  in 

Since  j  does  not  vary  linearly  with  the  uncertainty  at  high  and  low 

was  considered.  This  uncertainty  was  determined  by  utilizing  Figure  19 
and  the  uncertainty  in  Maximum  Slope  as  2.0%. 

=  3.5.  X=  0  =  +  6. OX 

tu 

N  =  30,  X=  .04  ‘^^u  =  +  9.6% 

tu  ’  rr- —  — 


Returning  to  Equation  (1)  and  substituting  as  appropriate: 


For  high  N  ;  'tL.  ^  \ 

^  J  1 

|^(.ofe/+  (01)’-+ 01) V  ('-^Y 

]\  .0t3 

tVj'  r 

J 

For  low  N_:  =  I 

^  J  1 

t  4  (.  oi)’-  +(.  oo"^  + 

'4 

=  .o<?8 

J 

Following  the  above  type  of  analysis,  the  following  additional 
uncertainties  were  calculated, 

m  +1.0% 

Np  ±  2 . 5% 

f  +4.3% 


21 


8.  DISCUSSION  OF  RESULTS. 


The  60  X  60  mesh  screen  matrix  was  used  solely  as  a  test  matrix,  a 
matrix  whose  heat  transfer  and  friction  characteristics  were  well  docu¬ 
mented.  Figure  7  and  Table  I  present  the  results.  There  is  quite  good 
agreement  in  the  friction  factor  but  a  slight  variation  in  j  factor.  This 
difference  is  attributed  to  the  difference  in  the  geometries  of  the  test 
matrices.  Exact  reproduction  of  a  screen  matrix  is  difficult  to  achieve. 
The  correlation  was  such  that  the  experimental  test  setup  was  considered 
capable  of  delivering  reliable  data. 

For  further  verification  on  the  reliability  of  the  test  setup.  Figures 
8  to  10  and  Tables  II  to  IV  present  the  test  results  for  the  Plate-Fin 
Reference  Matrix,  the  10°  Skew  Stainless  Steel  Matrix  and  the  20°  Skew 
Stainless  Steel  Matrix  used  by  Howard  (6).  Although  this  was  not  the 
original  purpose  for  utilizing  these  matrices,  the  results  do  indicate 
conclusively  that  the  "integral  heater"  is  as  good  or  better  than  the 
"sliding  drawer"  technique. 

The  original  purpose  for  utilizing  these  matrices  was  to  try  and  ex¬ 
tend  Howard's  curves  into  the  higher  Reynold's  number  region  (lower 
region).  The  points  obtained  in  this  region  are  plotted  on  the  curves  and 
the  perinent  data  are  recorded  in  the  Tables.  None  of  them  follow  the 
straight  line  of  their  predecessors.  There  are  at  least  two  reasons  for 
this  obvious  departure  from  the  norm.  First,  it  is  a  well  known  fact 
that  the  "Maximum  Slope"  Method  has  an  inherently  large  error  in  the  low 
N^^  region  of  1.5^^^<3.5.  This  fact  was  originally  published  by  Locke 
(1)  for  0.  Howard  (5)  took  into  account  longitudinal  conduction  and 
published  an  error  curve  that  includes  the  effect  of  longitudinal  conduc¬ 
tion  on  an  error  in  see  Figure  19.  In  the  regions  below  <1.5 
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the  error  is  acceptable,  yet,  a  discrepancy  still  exists.  This  points 
out  the  second  reason  for  the  discrepancy.  The  maximuni  slope  at  high 
flow  rates  (low  occurs  at  time  t  =  0+,  Instrumentation  prevents 

accurate  measurement  of  the  true  maximum  slope  at  this  time  due  to  in¬ 
herent  delay  or  lag  time. 

Extrapolation  of  the  time- temperature  curves  back  to  time  t  =  0  was 
attempted  in  an  effort  to  remove  the  effects  of  instrumentation.  The 
method  used  was  either  incorrect  or  unreliable.  Satisfactory  results 
were  unobtainable. 

The  merits  of  the  plate-fin,  10®  skew  and  20°  skew  matrices  as  heat 
transfer  surfaces  was  discussed  thoroughly  by  Howard  in  reference  (6). 

Figures  11,  13  and  14  and  Tables  V,  VI  and  VII  present  the  results 
for  a  new  geometry  tested  and  a  new  idea  in  material.  Figure  11  and 
Table  V  present  the  results  for  a  20°  skew  nickel  matrix  constructed  from 
slotted  plates.  Appendix  C  explains  fully  the  geometry  of  the  slots  and 
the  calculations  that  determine  the  basic  geometric  parameters  p,  ^  , 
and  r^.  The  effects  of  longitudinal  conduction  was  considered  negligible 
as  the  matrix  was  constructed  so  that  the  air  flow  was  normal  to  the  long 
dimension  of  the  slots.  Plotting  all  the  points  uncorrected  for  conduc¬ 
tion  showed  an  erratic  departure  from  the  expected  results  for  200. 

It  was  then  decided  that  perhaps  the  basic  assumption  was  not  warranted 
completely.  Corrections  for  longitudinal  conduction  were  made  where 
significant  and  these  points  were  plotted.  These  corrected  values  of  j 
are  overestimations,  for,  although  is  corrected  for  the  effects  of  the 
slots,  the  matrix  length  L  is  not  the  true  conduction  path  length  but  at 
best  a  lesser  approximation.  Since  X  is  obviously  not  equal  to  zero, 
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the  true  value  of  j  is  probably  between  the  uncorrected  j  (  A =  0)  and 
the  j  corrected  for  longitudinal  conduction  using  L  as  an  approximation 
for  conduction  path  length.  This  was  considered  in  the  interpretation  of 
results . 

Figure  12  compares  the  20°  skew  stainless  steel  matrix  with  the  per¬ 
forated  nickel  matrix.  There  is  a  70  to  1007,  increase  in  j,  while  f  in¬ 
creases  10  to  257,. 

It  must  be  understood  that  these  two  matrices  are  geometrically  similar 
only  in  overall  dimensions.  Actually,  the  heat  transfer  areas  were  cal¬ 
culated  on  different  bases.  The  20°  skew  stainless  steel  heat  transfer 
area  was  based  on  the  total  surface  area  available  for  heat  transfer  while 
the  heat  transfer  area  of  the  20°  skew  perforated  was  based  on  this  same 
area  less  the  area  of  the  slots.  Also,  not  included  in  the  heat  trans¬ 
fer  area  is  the  peripheral  area  of  the  slots  normal  to  the  flow.  If 
the  20°  skew  perforated  j  factor  were  based  on  the  total  available  area 
for  heat  transfer,  the  values  plotted  would  be  decreased  approximately 
40%.  The  perforated  material  still  represents  a  substantial  improvement 
in  the  performance.  This  improvement  is  attributed  to  the  reduction  of 
inefficient  stagnant  areas  within  the  matrix  and  the  alteration  of  the 
thermal  and  hydrodynamic  boundary  layers  caused  by  the  slots. 

Figure  13  and  Table  VI  present  the  results  for  a  "parallel  plate"  mat¬ 
rix  constructed  from  brass  shim  stock.  As  can  be  seen  in  Figure  13,  the 
heat  transfer  performance  falls  far  below  that  of  theoretical  laminar 
flow  performance  for  parallel  plates.  Several  things  may  account  for  this: 
the  aspect  ratio  of  the  flow  channels  in  this  matrix  was  77  vice  co  for 
truly  parallel  plates,  and  also,  the  flow  channels  were  not  uniform  in 
cross-section  from  one  to  another.  The  non-uniformity  in  cross-section 
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results  in  uneven  mass  flow  through  the  matrix,  hence  inefficient  use 
of  the  available  heat  transfer  surfaces,  and  therefore  a  decrease  in  j. 

The  decrease  in  friction  factor  is  possibly  due  to  the  sharp  angle 
regions  where  the  formed  sections  contact  the  flat  plate.  The  wall  shear 
stress  in  this  type  of  region  is  small  which  would  tend  to  reduce  the 
value  of  f.  These  regions  would  also  tend  to  produce  stagnant  regions 
previously  mentioned,  having  a  diminishing  effect  also  on  j. 

Figure  14  and  Table  VII  present  the  results  for  a  "parallel  plate" 
matrix  constructed  from  perforated  nickel  plate.  This  matrix  was  construct¬ 
ed  in  the  identical  fashion  as  the  brass  matrix.  Although  it  is  not  read¬ 
ily  apparent  in  Figure  14,  the  heat  transfer  characteristics  have  been 
improved  nearly  807o,  while  the  friction  factor  has  remained  the  same. 

Again  it  is  necessary  to  m.ake  note  of  the  fact  that  the  heat  transfer  area 
of  the  two  matrices  were  calculated  on  different  bases.  These  bases  are 
the  same  as  discussed  for  the  20°  skew  matrices.  Applicable  also  is  the 
approximate  reduction  by  40%  of  the  plotted  values  of  the  j  factor  if 
all  of  the  area,  in  the  perforated  material,  available  for  heat  transfer, 
were  considered  in  the  calculations. 

Figure  15  is  a  comparison  of  the  two  "parallel  plate"  matrices.  The 
improvement  realized  by  using  perforated  material  is  readily  apparent, 
even  considering  the  40%  reduction  mentioned  above.  The  additional  per¬ 
formance  gain  in  the  transferral  of  heat  is  again  attributed  to  the  re¬ 
lieving  of  the  stagnant  areas  and  the  alteration  of  the  boundary  layers 
caused  by  the  slots. 

The  fact  that  the  friction  factor  did  not  increase  does  not  lend 
itself  to  previous  conclusions  (20°  skew) ;  that  the  presence  of  slots 
increases  the  friction  facto) ,  Some  additional  factor  must  be  necessary. 
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’T’erhapp,  the  skewing  coupled  with  the  slots  contribute  more  to  friction 
than  straight,  slotted  channels. 

In  the  determination  of  the  j  factor  for  the  perforated  "parallel 
place"  matrix  the  conduction  parameter,  A  ,  was  considered  zero.  Here 
again  the  plotted  points  departed  from  the  expected  straight  line  indicat¬ 
ing  that  X  ,  in  fact,  was  not  zero.  Corrections  for  longitudinal  conduc¬ 
tion  were  made  using  the  matrix  length,  L,  as  an  approximation  for  the 
conduction  path  length.  These  points  were  plotted  also.  Since  the  cor¬ 
rected  points  are  overestimations  and  A  /  0,  the  true  value  of  j  lies 
between  these  limits. 

The  designer  of  a  compact  heat  exchanger  is  generally  restricted  as  to 
frontal  area,  volume,  i^eight  and  pressure  drop.  There  is  presently  coming 
into  use,  several  means  of  comparing  heat  transfer  surfaces  with  these 
restrictions  in  mind.  Essentially,  various  parameters  are  compared 
graphically  so  that  the  relative  merit  of  the  various  surfaces  may  be  es¬ 
tablished. 

The  ratio  j/f  is  plotted  vs.  Reynold's  Number  as  a  means  of  relating 
frontal  flow  area  with  pressure  drop.  For  a  given  pressure  drop  a  high 
j/f  characteristic  implies  a  small  flow  frontal  area,  j/f  when  plotted 
vs.  is  known  as  the  Flow  Frontal  Area  "Goodness"  factor. 

The  heat  transfer  power,  is  plotted  vs.  the  flow  friction  power, 

^STD’  resultant  characteristic  curve  is  a  measure  of  relative  core 

volumes.  Generally,  the  surface  with  a  high  h  vs.  E  characteristic 

b  1  D  b  1  D 

will  tend  to  have  a  small  core  volume,  plotted  in  this  fashion  is 

known  as  the  Volume  "Goodness"  factor. 

Figure  16  compares  the  20“  skew  matrices  and  the  "parallel  plate" 
matrices.  Although  the  "Goodness"  factors  are  not  infallible,  they  are 
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Indicative  of  relative  performance.  Tn  this  light,  the  perforated 
material  should  be  excellent  matrix  material. 
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9. 


CONCLUSIONS. 


Based  upon  the  correlation  of  present  data  with  that  obtained  by 
previous  experimenters,  the  experimental  test  rig  design  presented  here¬ 
in  is  considered  capable  of  producing  reliable  test  data  on  the  heat  trans¬ 
fer  and  friction  characteristics  of  a  core  module. 

Basic  heat  transfer  and  flow  friction  characteristics  for  a  new 
geometry  and  a  new  idea  in  material  are  presented  in  Figures  11  to  15. 

The  geometric  and  physical  properties  are  summarized  in  Figures  5  and  6. 
There  appears  to  be  a  definite  advantage  to  constructing  matrices  from 
perforated  material.  Performance  is  improved  considerably  and  based  on 
the  "Goodness”  factors  (Figure  16)  size  and  weight  advantages  are  realiz¬ 
ed. 

The  20°  skew  perforated  nickel  matrix  was  the  top  surface  in  perfor¬ 
mance.  This  matrix  had  the  highest  Flow  Frontal  Area  "Goodness"  factor, 
as  well  as  the  highest  Volume  to  surface  "Gcodness"  factor.  This  does  not 
mean  that  this  is  the  "optimum"  surface  for  any  given  application  but  it 
does  merit  consideration. 

In  both  geometries  tested  (20°  skew  and  "parallel  plate")  the  per¬ 
forated  material  out-performed  its  counterpart.  This  would  also  Indicate 
that  further  testing  of  this  type  of  material  (slotted)  is  warranted. 

It  would  seem  that  this  material  and  some  geometry  should  result  in  an  out¬ 
standing  heat  exchanger  surface. 
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10,  RECOMMENDATIONS  FOR  FUTURE  WORK. 


1)  Some  means  should  be  found  to  reliably  determine  N^.^  in  the 

high  range  (1.0  ^  3.5),  Most  promising  is  the  use 

of  a  periodic  or  cycling  technique  in  which  effectiveness  is 

measured  and  used  to  find  N^  rather  than  the  maximum  slope. 

tu  ^ 

This  ability  is  necessary  if  performance  curves  are  to  be  extend¬ 
ed  to  500. 

2)  More  work  can  be  done  with  perforated  material.  Different 
geometries  should  be  tested  to  see  if  consistently  high  per¬ 
formances  are  realized  as  they  have  been  in  the  past.  Addition¬ 
ally,  this  type  of  surface  is  suspect  of  fouling.  If  these 
materials  are  as  good  as  they  seem  to  be,  then  examination  of 
the  effects  of  fouling  on  their  performance  is  in  order. 

3)  As  more  surfaces  become  available,  they  should  be  tested 

and  their  heat  transfer  and  friction  characteristics  determined. 
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FIGURE  1.  SCHEMATIC  DIAGRAM  OF  TEST  APPAmTUS 


18-8  St. 

Matrix  Material  Steel 

Specific  Heat  (c^)  BTTJ/lb  »F  0-12 

Thermal  Ocnductivity  (kg)  BTU/hr  -P  ft  10.0 


Material  Thickness,  inches  ^ 

Total  Heat  Transfer  Area  (A)  ft 
Frontal  Area  (A^^)  ft  ^ 

Total  Conduction  Area  (A^)  ft  ^ 
Plow  Cross  Section  Area  (A^)  ft 
Matrix  Volume  (V^)  ft  ^ 
Matrix  Density  Ih/ft 


0.017 

4.0159 

0.07055 

0.02215 

0.0484 

0.00199 

155. 


Hydraulic  Diameter  (D^^)  ft 
Compactness  (y0)  ft^/ft^ 
Porosity  (p) 


0.00156 

2009. 

0.686 


PIGTJRH  3 


Physical  and  Geometric  Properties  of 
60x60  Mesh  Screen  Matrix 
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Type  302 


Matrix  Material 

S.  Steel 

Specific  Heat  (c^)  BTU/lh  op 

0.11 

Thermal  Conductivity  (k^)  BTU/hr  op  ft 

7.0 

Material  Thickness,  inches 

0.001 

Total  Heat  Transfer  Area  (A)  ft^ 

50.12 

Frontal  Area  (A^^)  ft^ 

0.07023 

Total  Conduction  Area  (A^)  ft^ 

0.00752 

Plow  Cross  Section  Area  (A^)  ft^ 

0.06271 

Matrix  Volume  (V^)  ft^ 

0.01170 

Matrix  Density  (p^)  Ib/ft^ 

52.9 

Hydraulic  Diameter  (Dj^)  ft 

.001388 

Compactness 

2573. 

Porosity  (p) 

0.893 

FIGURE  4  Physical  and  Geometric  Properties  of 
Plate-Pin  Reference  Matrix 


35 


Matrix  Material 

Specific  Heat  (dg)  BTU/lb 

Thermal  Cond.  (kg)  BTU/hr  ep  ft 

Material  Thickness,  inches 

2 

Total  Heat  Transfer  Area  (A)  ft 

Frontal  Area  (A-„)  ft^ 

fr  P 

Total  Conduction  Area  (A„)  ft 

B  2 

Flow  Cross  Section  Area  (A^)  ft 
Matrix  Volume  ft^ 

Matrix  Density  Ih/ft^ 

Hydraulic  Diameter  (Dg)  ft 
Compactness  (/3)  ft^ft^ 

Porosity  (p) 


IQo  SKEW  20°  SKEW  20°  SKEW 


Type  502 
S. Steel 

Type  502 
S. Steel 

vperf . 

Nickel 

0.11- 

0.11 

0.1065 

7.0 

7.0 

58.7 

0.001 

0.001 

0.0022 

18.58 

18.04 

15.41 

0.07024 

0.0702 

0.0754 

0.00464 

0.00452 

0.00549 

0.0656 

0.0657 

0.06412 

0.0117 

0.01157 

0.01211 

52.4 

52.22 

50.1 

.002581 

.002416 

.002528 

1570. 

1549. 

1272. 

0.954 

0.956 

0.875 

PIGUEE  5  Physical  and  Geometric  Properties  of 
Skew  Matrices 
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Matrix  Material 

70-30 

Brass 

Perf. 

Nickel 

Specific  Heat  (Cg)  BTU/lb  op 

0.092 

0.1065 

Thermal  Conductivity  Ckg)  BTU/hr  op  ft 

66. 

38.7 

Material  Thickness,  inches 

0.001 

0.0022 

Total  Heat  Transfer  Area  (A)  ft^ 

22.057 

18.586 

Frontal  Area  (A^^)  ft^ 

0.06986 

0.0695 

Total  Conduction  Area  (Ag)  ft^ 

0.00551 

0.00416 

Plow  Cross  Section  Area  (A^)  ft^ 

0.06434 

0.0584 

Matrix  Volume  (V^)  ft^ 

0.01164 

0.01158 

Matrix  Density  Ih/ft^ 

41.5 

63.6 

Hydraulic  Diameter  (D^^)  ft 

.001944 

.001928 

Compactness  (^)  ft^/ft^ 

1895. 

1605. 

Porosity  (p) 

.921 

.840 

FIGURE  6  Physical  and  Geometric  Properties  of 
"Parallel  Plate"  Matrices 
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39 


Figure  9  10°  Skew  Stainless  Steel  Matrix 

I 
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Figure  10  20®  Skew  Stainless  Steel  Matrix 

I 
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Figure  12  Comparison  of  perforated 

and  non-perforated  20°  Skew 
.  Matrices 
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Figure  15  ."Parallel  Plate"  Brass  Llatrix 
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Figure  14  "Parallel  Plate"  Nickel  (perforateji)  Matrix 

)  .  • 
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Pigure  15  Comparison  of  Perforated  and 

non-perf orated  "Parallel  Plate" 
I  .  Matrices 
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200 


0.1 


(T)20o  Skew  Perforated  Nickel 
@ "Parallel  Plate"  Perforated  Nickel 
@  20°  Skew  Stainless  Steel 
"Parallel  Plate"  Brass 


0.5 

0.3 


Pigure  16  Design  "Goodness"  Factors 
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TIME 

FIGURE  17  GENERALIZED  COOLING  CURVE 

NOTE:  All  measurements  are  made  In  inches 
to  facilitate  computations. 
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NUMBER  OF  TRANSFER  UNITS  VERSUS  MAXIMUM  SLOPE 
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Figure  18  Numlaer  of  Transfer  Units  (N. 
Maximum  Slope 
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Maxiiirutm  Slope  vs 
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FIGURE  22,  Matrix  Holder  modified  with  probe 
for  velocity  profile  measurement 
and  temperature  distribution 
measurement  (Downstream  View) 


FIGURE  23.  Nichrome  Wire  Heater  (Upstream  View) 


FIGURE  2U,  Nichrome  Wire  Heater  (Downstream  View) 


PROGRAM  MAXSLO 
L  =  2 
K=  1 

30  READ  300»XML,AS,RH,P0R,AFL, AHT,  WM 
300  FORMAT(7F10.5) 

40  READ  400,EK,CK,SK,CM,N 
400  FORMAT(4F10.5, 12) 

M=1 

10  READ  100,DO,DELPO,GAM,ATMP,COR, UH.UFR 

100  FORMAT(7F10.5) 

READ  101,CS,DELH  , HS , SLO , TEMPO 

101  F0RMAT(5F10.5) 

EMDOT=359.0*DO*DO*SQRTF(DELPO»GAM*ATMP/29.92)*COR 

CP=0.24 

RNUMP=5.09*EMDOT/UFR 
RNUMH=(4.0*XML*EMD0T)/( AHT*UH) 

RNUMFR=RNUMH»UH/UFR 

CAPFL=EMDOT*CP/3600.0 

CAPS=WM«CM 

SLOMAX= (CAPS •SLO) / (CAPFL*CS1 
CONDPAR=(SK*AS)/ (XML*EMCOT*CP) 

B=CK+EK 

C=1 .0  +  P0R*»2 

G=EMD0T/{AFL*3600.0) 

DELP=DELH*. 03613 

PS=HS*. 03613 

PA=ATMP».4912 

P1=PA-PS 

P2=P1-DELP 

PM=(Pl+P2)/2. 

RH0M={PM»144.0) /{53.3»( TEMPO+460.0) ) 

FF1=( (64*4*RH0M*DELP*144.0)/IG*G} ) 

FF2=C*DELP/PM 


FF3=FFl-FF2-8 
FFR=FF3*RH/XML 
WRITE  OUTPUT  TAPE 
F0RMAT(9H  RUN  NO. 
WRITE  OUTPUT  TAPE 


4,4,M 
I2///  ) 

_  _ _ _ _  4, 1 ,EMDOT,RNUMP,RNUMH,CAPFL,CAPS 

1  F0RMAT(7H  EMDOT= F 1 0. 5 , 3X , 7H  RNUMP=E 20. 8 , 3X , 7H  RNUMH=F10. 
17H  CAPFL=F10.5,3X,6H  CAPS=F1 0.5///) 

WRITE  OUTPUT  TAPE  4. 2 , CONDPAR , S LOMA X 

2  F0RMAT(9H  CONDP AR=F 1 0 . 5 , 3X , 8H  S LOMA X=F 1 0. 5/ // ) 

WRITE  OUTPUT  TAPE  4 , 3 , G  ,  RNUMFR , FFR 

3  F0RMAT(3H  G=F 1 0 .5, 3X , 8H  RNUMFR»F10. 5, 3X. 5H  FFR=F10.5/) 
M=M  +  1 

IF(N-M)50, 10,10 
50  K=K^^1 

IF(L-K)20,30,30 
20  STOP 
END 
END 


Figure  27  Program  MaxSlo  Readout 
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TABLE  I  SUTJIv’lARY  OF  HEAT  TRANSFER  AND  FRICTIONAL  RESULTS 
60x60  MESH  STAINLESS  STEEL  SCREEN  MATRIX 


HEAT  TRANSFER  RFSTTT.t?^. 

R^  REYNOLD’S 
/  NUMBER 

1  431.57 

2  234.81 

5  159.26 

4  129.69 

5  106.78 

6  79.88 

7  50.69 

8  29.60 

9  15.43 

10  8.49 


CONDUCTION 

MAXIMUM 

No, 

PARAMETER 

SLOPE 

tu 

0.0 

.58333 

3.12 

0.0 

.68566 

4.96 

0.0 

.74639 

6.09 

0.0 

.80395 

7.25 

0.0 

.82926 

7.77 

0.0 

.87817 

8.83 

0.0 

1.05579 

11.90 

0.0 

1.22105 

17.94 

0.0 

1.57029 

30.19 

0.0 

2.11897 

55.65 

.0267 

.0425 

.0521 

.0620 

.0659 

.0755 

.1040 

.1535 

.2584 

.4760 


ISOTHERMAL  FRICTION 


RUN  REYNOLD’S 

4  NUMBER 

1  437.93 

2  238.27 

3  161.60 

.4  131.60 

5  108.36 

6  81.07 

7  51.44 

8  30.04 

9  15.66 

10  8.62 


RESULTS ! 

FRICTION 

FACTOR 

.68548 

.82196 

.89263 

.97400 

1.06410 

1.23243 

1.57979 

2.31547 

4.47090 

9.50629 
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TABLE  II  SUraiAEY  OP  HEAT  TRANSFER  AND  FRICTIONAL  RESULTS; 
PLATE-PIN  REFERENCE  MTRIX 


HEAT 

TRANSFER 

RESULTS;. 

RUN 

/ 

REYNOLD’S 

NUTvIBER' 

% 

N  N 

jmst^pr 

0 

CONDUCTION 

PARAMETER 

X 

MAXIMUM 

SLOPE 

^tu 

•^tu 

(A=o) 

1 

505.17 

.00277 

.00128 

.52585 

1.68 

1.68' 

2 

433.26 

.00429 

,00150 

.55712 

2.60 

2.60- 

3 

387.19 

.00505 

.00167 

.58044 

5.06 

5.06 

4 

522.08 

.00654 

.00201 

.62585 

5.84 

•  3.84 

5 

248.92 

.00842 

.00260 

.69577 

5.10 

5.10 

6 

187.25 

.01162 

.00546 

.79434 

7.04 

7.04 

7 

125.74 

.01761 

.00515 

.94678 

10.67 

10.44 

8 

75.04 

.02889 

.00862 

1.16341 

17.51 

16.21 

9 

42.68 

.05699 

.01516 

1.43365 

34.54 

25.07 

10 

25.81 

.08405 

.02717 

1.43851 

50.93 

25.22 

ISOTHERJJAL  FRICTION  RESULTS; 


RUN 

/ 

REYNOLD’S 

NUMBER 

% 

FRICTION 

FACTOR 

f 

1 

512.59 

.02197 

2 

439.62 

.02485 

3 

392.86 

.02795 

4 

326.79 

.05304 

5 

252.57 

.04210 

6 

189.99 

.05433 

7 

127.58 

.08205 

8 

76.14 

.12974 

9 

43.30 

.25164 

10 

24.17 

.39657 
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TABLE  III  SmOIARY  OP  HEAT  TRANSFER  AND  FRICTIONAL  RESULTS: 
10°  SKEW  STAINLESS  STEEL  MATRIX 


HEAT 

TRANSFER 

RESULTS: 

RUN 

REYNOLD’S 

N  N 
"srTR 

CONDUCTION 

MAXI&IUM 

^tu 

# 

NUIffiER 

PARAMETER 

SLOPE 

% 

A 

(A=0) 

1 

805.27 

.00575  ‘ 

.00081 

.46555 

1.52 

1.52 

2 

645.52 

.00440 

.00101 

.51088 

1.55 

1.55 

5 

517.86 

.00695 

.00126 

^55129 

2.45 

2.45., 

4 

402.58 

.01075 

.00162 

.62114 

5.78 

5.78 

5 

505.00 

.01488 

.00214 

.70167 

5.24 

5.24 

6 

200,72 

.02271 

.00524 

.85894 

8.00 

7.98 

7 

119.58 

.05642 

.00544 

1.02840 

12.85 

12.46 

8 

68.74 

.06576 

.00955 

1.27448 

22.46 

19.62 

9 

57.61 

.10788 

.01728 

1.44695 

58.00 

25.55 

ISOTHERMAL  FRICTION  FACTOR: 


RUN 

/ 

REYNOLD’S 

NUIffiER 

Nr 

FRICTION 

FACTOR 

f 

1 

817.04 

.02791 

2 

652.71 

.05405 

5 

525.45 

.04092 

4 

408.46 

.04951 

5 

509.46 

.06545 

6 

205.79 

.09412 

7 

121.41 

.15158 

8 

70.00 

.27472 

9 

58.16 

.47276 
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TABLE  IV  SmmiARY  OP  HEAT  TRANSFER  AND  FRICTIONAL  RESULTS: 
20°  SKEW  STAINLESS  STEEL  MATRIX 


HEAT 

TRANSFER 

RESULTS:  , 

RUN 

# 

REYNOLD’S 

NUTffiER 

^R 

N  N 
^^sr  PR 

CONDUCTION 

PARAMETER 

A 

MAXIMUM 

SLOPE 

^tu 

(A=0) 

1 

855.38 

.00594 

.00078 

.47559 

1.56 

1.56 

2 

728.65 

.00414 

.00090 

.48768 

1.45 

1.45 

3 

636.57 

.00475 

.00105 

.52594 

1.64 

1.64 

4 

535.10 

.00805 

.00122 

.56474 

2.78 

2.78 

5 

412.85 

.01182 

.00158 

.65657 

4.09 

4.09 

6 

512.06 

.01524 

.00209 

.70256 

5.26 

5.26 

T 

208.65 

.02167 

.00515 

.81581 

7.48 

7.48 

8 

126.47 

.05676 

.00517 

1.02539 

12.69 

12.55 

9 

71.21 

.06593 

.00917 

1.28940 

22.76 

20.09 

10 

59.79 

.09441 

.01646 

1.39790 

52.59 

25.79 

ISOTHERMAL  FRICTION  FACTOR: 


RUN 

/ 

REYNOLD’S 

NUMBER 

^R 

FRICTION. 

FACTOR 

f 

1 

848.11 

.05197 

2 

739.39 

.05471 

3 

645.95 

.05908 

4 

542.97 

.04465 

5 

418,94 

.05451 

6 

316.65 

.06856 

7 

211.73 

.10255 

8 

128.53 

.16015 

9 

72.26 

.28784 

10 

40.40 

.  .48428 

table  V  SUM,IARY  OP  HEAT  TRANSPER  AITO  PRICTIOHAL  RESULTS 
20®  SKEW  NICKEL  MATRIX  (PERFORATED  PLATES) 


heat  transfer  RESm.TS. 


RUN  REYNOLD'S 
^  NUI/CBER 
% 

1  966.07 

2  835.43 

5  738.42 

4  619.79 

5  484.37 

6  362.97 

7  239.79 

8  145.03 

9  82.26 

10  46,22 

11  23.50 


N  N  ^ 

sr^R 

CONDUCTION  MAXDRJM 

N^ 

N 

PARAJTETER 

X 

SLOPE 

tu 

tu 

(A=cO) 

.00900 

.00337 

.56394 

2.72 

2.72 

.01146 

.00390 

.60157 

3.42 

3.46 

.01219 

.00441 

,  .61776  • 

3.70 

3.68 

.01418 

.00526 

.64762 

4.24 

4.28 

.01802 

.00674 

.71212 

5.41 

5.44 

.02497 

.00898 

.81013 

7.58 

7.54 

.03517 

.01359 

.95653 

11.29 

10.62 

.05199 

.02243 

1.14669 

19.65 

15.70 

.a7879 

.03951 

1.39808 

23.79 

.08147 

.07039 

1.42079 

24.60 

.09134 

.13843 

1.50225 

27.58 

ISOTHER]\IAL  FRICTION  ffAnrpm?. 


RUN 

REYNOLD'S 

FRICTION 

/ 

NUIffiER 

FACTOR 

% 

f 

1 

980.29 

.03316 

2 

847.72 

.03709 

3 

749.20 

.04116 

4 

628.91 

.04686 

5 

491.63 

.05802 

6 

368.31 

.07208 

7 

243.32 

.10258 

8 

147.16 

.16737 

9 

83.46 

.29667 

10 

46.89 

.49889 

11 

23.85 

1.09090 
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TABLE  VI  SUMMARY  OF  HEAT  TRANSFER  AND  FRICTIONAL  RESULTS’: 
'•PARALLEL  PLATE"  BRASS  MATRIX 


HEAT  TRANSFER  RESULTS: 


RUN 

/ 

REYNOLD’S 

NUMBER 

% 

N  N 
"sr  PR 

CONDUCTION 

PARAMETER 

X 

BiAXIMUM 

SLOPE 

^tu 

(X=0) 

1 

683.97 

.00376 

.00896 

.49410 

1.44 

1.46 

2 

551.23 

.00527 

.01112 

.54861 

2.27 

2.40 

5 

437.55 

.00748 

.01401 

.59598 

3.22 

3.36 

4 

340.15 

.01038 

.01802 

.66717 

4.47 

4.63 

5 

255.77 

.01335 

.02390 

.73472 

5.75 

5.63 

6 

168.67 

.02148 

.03634 

.87166 

9.25 

8.71 

7 

102.18 

.03483 

.05999 

.97562 

15.00 

11.15 

8 

57.06 

.  - 

.10761 

1.07854 

- 

- 

'9 

33.37 

- 

.18369 

.98829 

- 

ISOTHERMAL  FRICTION  FACTOR; 


RUN 

# 

REYNOLD’S 

NUrfflER 

^R 

FRICTION 

FACTOR 

f 

1 

694.19 

.02846 

2 

559.46 

.03454 

3 

444.09 

.04292 

4 

345.24 

.05240 

5 

259.52 

.06907 

6 

171.19 

.10680 

7 

103.71 

.17239 

8 

57.90 

.32370 

9 

33.86 

.53085 
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TABLE  VII  SmC-lAHY  OP  HEAT  TRANSFER  AND  FRICTIONAL  RESULTS: 

"PARALLEL  PLATE"  NICKEL  MATRIX  (PERFORATED  PLATES) 


HEAT 

TRANSFER 

RESULTS: 

RUN 

# 

REYNOLD* 

NUMBER 

^R 

S  N  N 

d  . 

CONDUCTION 

PARAMETER 

A 

MAXIMUM 

SLOPE 

^tu 

^tu 

(A=o) 

1 

819.49 

.00590 

.00393 

.54781 

2.32 

2.36 

2 

647.21 

.00810 

.00497 

: 58864 

3.18 

3.24 

3 

518.89 

.00988 

.00620 

.62923 

3.89 

3.93 

4 

400.50 

.01458 

.00804 

.73180 

5.77 

5.83 

5 

304.49 

.01891 

.01058 

.81941 

7.61 

7.56 

6 

200.68 

.02659 

.01602 

.95459 

11.31 

10.63 

7 

121.15 

.04039 

.02658 

1.16142 

21.61 

16.15 

8 

67.53 

.05502 

.04768 

1.34786 

- 

22.00 

9 

37.92 

.04945 

.08478 

1.27951 

- 

19.70 

ISOTHERMAL  FRICTION  FACTOR: 


RUN 

# 

REYNOLD’S 

NUT/BER 

FRICTION 

FACTOR 

f 

1 

831.56 

.02395 

2 

656.72 

.02955 

3 

526.51 

.03577 

4 

406.49 

.04439 

5 

309.04 

.05682 

6 

203.63 

.08725 

7 

122.96 

.14242 

8 

68.54 

.26776 

9 

38.48 

.4587f 
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APPENDIX  A. 


DATA  REDUCTION  RELATIONS 


The  following  summary  of  some  of  the  data  reduction  relations  is 
given  to  clarify  some  of  the  definitions  of  the  geometrical  parameters, 
to  show  how  they  influence  the  deduced  results  and  to  facilitate  an 
understand  as  to  the  exact  procedure  required  to  produce  from  raw  data 


a  point  on  the  j  vs  N^  and  f  vs  N^^  curves 


The  three  most  significant  geometric  parameters  that  afford  a  means 
of  comparison  between  matrices  are  so  defined  that  there  exists  a  unique 
relation  between  them.  Independent  measurement  of  any  two  establishes 
the  third  or  serves  as  a  check  on  its  independent  measurement.  The 
three  primary  surface  geometry  parameters  of  a  matrix  are: 

1)  HYDRAULIC  DIAMETER: 


(A-1) 


2)  POROSITY: 


Matrix  Flow 
Void  Volume 


A 


c 


p  =  Total  Volume 


(A-2) 


A. 


‘fr 


3)  AREA  COMPACTNESS:  Total  Heat  Transfer 

5iiirfano 


(A- 3) 


NOTE:  If  equation  (A-2)  is  divided  by  (A-3)  the  unique  relation  mentioned 
previously  is  formed. 


(A-4) 


MAXIMUM  SLOPE; 

As  previously  discussed,  through  the  work  of  Locke,  Mondt,  Crestwick 
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and  .oward,  it  has  been  established  '.h;' :  the  axl  slope  jf  a  general¬ 


ized  heating  or  cooling  curve  is  a  unique  function  of  and  /A 
conduction  parameter, 


*'f2  ~  ^i 

fi  i 


=  X  ) 


T  =  generalized  time  variable 


T 


:  air  temp,  downstream 


T/n  , 


me  /W  c  © 
p  s  s 


and  d(  T  /N^^)  =  mc^/  (d  ©  ) 


f2  fl  +  1 
’^fl  ■  ^ 


t  „  -  t, 

and;  d(-^ - 

fl  ■  i 

combining  (A-5)  and  (A-6) 


t,,). 


d(  ■ 


fl 


d(  T/N,,  ) 


d(  ■ 


,  the 


(A-5) 


(A-6) 


(A- 7) 
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The  matrix  Is  heated  to  a  uniform  temperature,  t^.  The  heat  source 
is  secured  and  the  difference  between  downstream  and  upstream  temperatures 
(tj2  ■  recorded  vs.  time.  The  maximum  rate  of  change  of  this 

difference  with  time  may  be  determined  from  the  record. 

From  Figure  17 ; 


1 

X 


_ X 

Chart  speed 


dO  sec 


y  =  d(tf2  -  tfj) 


inches 


M  =  (t^j^  -  t^)  inches 


W  c  Matrix  heat  Capacitance 

-  =  Fluid  heat  Capacitance 

mCp 


TT  BTU/°F 
^  BTU/^F  sec 


Finally,  Equation  (A-7)  becomes: 


d(  ^f2  '  ^i  ) 
^fl  '  *^i 

d(  T 


max 


_L 

M 


Chart  speed 


With  this  value  of  the  maximum  slope  and  a  calculated  value  of  A 
Howard's  curves  or  tables  may  be  entered  to  obtain  a  value  of  N 

tu 


MASS  FLOW  RATE: 

The  mass  flow  rate  is  calculated  as  outlined  in  reference  (7): 
m  =  359  K  Dq  FaY  |  A  t 
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where 


K  =  flow  coefficient,  velocity  of  approach  factor 


Fa  =  thermal  expansion  factor 

Y  =  expansion  factor 

D  =  orifice  diameter,  inches 
o  5 

=1  specific  weight  of  fluid  flowing  at  the  orifice 

3 

temperature  lb/ ft  ,  corrected  to  existing 
atmospheric  pressure 

A =  Pressure  drop  across  theorifice  ,  inches  H2O 

REYNOLD’S  NUMBER: 


By  definition,  the  Reynolds  number  is 


R  JJL 

T'^here  G  is  the  mass  flow  velocity  based  on  the  flow  mean  cross- 
sectional  area,  A  . 


(A-8) 


Note,  if  (A-9)  is  substituted  into  (A-8) 


(A- 9) 


N 


R 


4m 


( 


J5l 

f 


) 


(A- 10) 


but  from  (A-4)  ^ 

Therefore,  for  a  given  range  of  mass  flows,  the  Reynolds  number  range 
for  a  given  matrix  is  proportional  to  the  inverse  of  its  compactness. 


For  ease  of  calculation,  recall  from  (A-3) 


substituting  into  (A-11)  yields 


4Lm 

A/^ 


(A-12) 


N 


R 


THE  FANNING  FRICTION  FACTOR; 

The  following  sketch  describes  the  flow  system  to  be  considered. 
1 _ a _ t _ 2 


1 

1 

;  ^ 

7/7777777777777777777, 

\ 

I 

1  ^ 

1  ^ 

7777777777777777777777, 

1 

1 

1 

TTT/////////////////// 

i 

For  gas  flow  heat  exchanger  application,  the  pressure  changes 
from  section  1  to  a  and  from  b  to  2  are  very  small,  relative  to  the 
total  pressure,  and  therefore  and  ^ 

The  relation  for  the  flow  stream  pressure  drop  calculation  for 
most  heat  exchanger  cores  is  (8); 


AP  -  V7 

P.  P. 


entrance  flow  core 

effect  acceleration  friction 


exit 

effect 


r-  Y-  po  ^  p^) 


(A-13) 


Solving  (A-13)  for  f  and  recalling  from  (A-1)  A/A^  =  L/r^  also 

J. 

'V  =  ^  ,  we  find; 


(A-14) 


Now  the  perfect  gas  law  may  be  applied  since  testing  is  done  with  air  at 
moderate  temperatures  and  pressures.  Equation  (A-14)  becomes  for  the 
isothermal  case: 
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f- 


T  (A-15) 


p  =  and  T,  = 

/  n’ 


U  _L  RT 

where  ^  =  p  ^  -  2  ^1  "  "2 

The  first  term  is  by  far  the  largest  contributor  to  the  friction  factor 


for  small  pressure  differentials,  therefore  making  the  approximations 


Equation(A-15)  becomes 


i-- 


where : 


(A-16) 


Kc  =  entrance  coefficient 
Ke  =  exit  coefficient 

whose  values  are  determined  from  appropriate  curves  in  Chapter  IV,  refer¬ 
ence  (8). 

This  is  the  equation  used  to  determine  the  friction  factor  as 
presented  in  this  thesis. 

If  a  further  assumption  is  made  such  that 

■f  »  ^  (A-U) 


and  equations  (A-4)  and  (A-9)  are  substituted.  Equation  (A-17)  becomes; 


•  ^  ^ ^  P  /  Af.r  \  p>^ 

V  L  /  /3 


(A-18) 
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Therefore,  for  a  given  matrix  the  friction  factor  is  directly  pro¬ 
portional  to  the  cube  of  the  porosity  and  indirectly  proportional  to  the 
compactness . 

f  (A- 19) 


COLBURN  i  FACTOR; 

We  have  A  and  maximum  slope  with  which  we  get  N^^  from  Table  VIII 
or  Figure  18  taken  from  reference  (5). 

The  j  factor  is  defined  as: 


j  =  2/3 

ST  PR 


G  c 


(A- 20) 


substituting  (A-9)  for  G  and  multiplying  by  /A 


1  =  N 

^  mr  A  P 


2/3 


(A-21) 


by  definition; 


hA 

ifiCp 


Ar  2/T 

Therefore:  j  =  —  Np^  (A-22) 

Note,  from  (A-1)  Ac/A  =  ’^H/L,  substituting  f  f ^  for 

r^  from  (A-4)  gives: 


2/3 


(A-23) 


which  leads  to  the  relation  that  for  a  given  matrix,  j  is  directly  pro¬ 
portional  to  the  porosity  and  inversely  proportional  to  the  compactness. 

J  OC  'V(S 

HEAT  TRANSFER  POWER  AND  FLOW  FRICTION  POWER; 

The  heat  transfer  power  per  unit  area  per  degree  temperature  differ¬ 
ence  is  derived  as  follows: 
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j  = 

ST  1 


2/3 


Solving  for  h: 


j  = 


k 


\ 

2/3 


p  M 


2/3 


Now  if  c^,  ^  and  are  evaluated  at  standard  properties  which  by 
convention  are  taken  as  those  for  dry  air  at  500 °F  and  one  atmosphere 

c  =  0.2477  BTU/lb  “F 
P 

/X  =  0.0678  Ib/hr  ft 


P 


=  0.0413  Ib/ft 


Thus,  the  heat  transfer  power  per  unit  area  per  degree  temperature 


hsTD  =  -02195  (  ^  ®TU/hr  ft^^F 

The  flow  friction  power  per  unit  area  is  derived  as  follows: 


Al.1 
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£  -^P(^)v(isS)  .  -f 


where  Ca  s  p  ^ 

E  =  -f 


4  (?:?)“  (fV*'] 


Now  if  ^  and  jO  are  evaluated  at  standard  properties  (previously 

defined)  the  flow  friction  power  per  unit  area  evaluated  at  the  above 

standard  conditions,  becomes: 

’  STD’ 

"sTD  =  ‘-ll  -  1°''  <  4^>  '  f 

Figure  16  provides  a  comparison  of  the  20®  skew  and  the  "parallel 
plate"  matrices  from  a  heat  transfer  power  vs.  friction  power  point  of 
view  vs.  .  For  comparison  purposes  the  surface  geometries 

have  been  adjusted  to  a  common  hydraulic  diameter  of  4r^  =  2  x  10  ^ft. 
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APPENDIX  B. 


DESCRIPTION  OF  EQUIPMENT 

To  obtain  heat  transfer  data,  a  transient  technique  was  employed. 
Briefly,  this  method  consists  of  heating  the  test  matrix,  to  a  uniform 
temperature  and  then  subjecting  the  matrix  to  a  step  change  in  fluid 
flow  temperature.  The  fluid  temperature  downstream  of  the  matrix  is 
recorded  versus  time. 

Friction  factor  was  calculated  from  pressure  drop  data  obtained 
from  static  pressure  taps  located  in  the  test  section  immediately  up¬ 
stream  and  downstream  of  the  matrix. 

The  equipment  necessary  to  complete  the  experiment  falls  into  one  of 
the  following  systems. 

1)  Fluid  source 

2)  Flow  metering  system 

3)  Temperature  measuring  system 

4)  Pressure  measuring  system 

5)  Fluid  heater 

6)  Matrix  holder  &  test  section  casing 

FLUID  SOURCE; 

The  working  fluid  was,  of  course,  air.  Air  was  provided  to  the 
apparatus  by  placing  the  entire  rig  on  the  inlet  to  a  lOOcfm,  6  stage 
centrifugal  air  compressor.  The  compressor  requires  an  operating  voltage 
of  120/  at  90amp. 

FLOW  METERING  SYSTEM;  (See  Figure  1) 

The  flow  metering  system  consists  of  an  orifice  in  a  3.08"  inside 
diameter  aluminum  tube.  The  orifice  section  was  constructed  in  accord¬ 
ance  with  the  ASME  specifications  of  reference  (7)  for  standard  D  and 
D/2  pressure  taps.  A  large  flow  range  was  obtained  by  means  of  7  orifice 
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plates  of  throat  diameter  2.310,  1.971,  1.232,  1.540,  .462  and  .308 
inches,  respectively.  Accurate  control  of  the  pressure  drop  across  the 
orifice  was  maintained  by  a  large  sliding  valve  downstream  of  the  ap¬ 
paratus  on  the  inlet  to  the  compressor  and  another  sliding  valve  on  the 
discharge  side  of  the  compressor.  There  is  an  additional  valve  on  the 
inlet  side  of  the  compressor  that  is  fully  open  or  fully  closed.  It  is 
normally  opened  at  flow  rates  less  than  150  Ib/hr  to  prevent  stalling  the 
compressor. 

TEMPERATURE  MEASURING  SYSTEM:  (See  Figure  1) 

Without  exception,  all  temperatures  were  measured  with  iron-constan- 
tan  thermocouples.  A  single  thermocouple  measured  the  orifice  tempera¬ 
ture  (T^)  just  upstream  of  the  orifice  at  the  inlet  to  the  apparatus. 

This  temperature  was  read  in  millivolts  on  a  RUBICON  CO.  portable  pre¬ 
cision  potentiometer,  converted  to  temperature  and  recorded  for  each  run. 
Measuring  is  a  group  of  5  thermocouples  in  series,  bound  together  and 
insulated  from  each  other  by  teflon  tape  and  inserted  in  a  1/8”  diameter 
aluminum  tube.  The  tube  was  placed  in  the  air  stream  with  a  small  window 
cut  in  it  facing  upstream.  At  this  point  ambient  air  temperature  was 
measured.  T2  and  T^  are  5  wire  thermocouple  grids  permanently  mounted  in 
the  test  section  cabinet.  Their  wires  are  connected  in  series  to  magnify 
the  thermocouple  output  voltage.  T^  is  a  movable  5  wire  thermocouple 
grid  located  in  the  matrix  holder  and  is  placed  adjacent  to  and  down¬ 
stream  of  the  matrix,  it’s  wires  are  also  in  series.  It’s  position  is 
apparent  in  Figure  21,  The  thermocouples  arranged  in  this  manner  provided 
flexibility  and  complete  coverage  of  desired  points  of  interest.  If  the 
thermocouple  output  voltages  from  T2  and  T^  are  bucked  against  each  other, 
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the  temperature  measuring  system  response  curves  may  be  generated  on  a 
Leeds  &  Northrup  "Speedomax"  temperature  recorder.  T2  vs  gives  the 
indication  of  temperature  uniformity  across  the  matrix  and  vs 
records  the  rate  of  change  of  the  difference  between  upstream  and  down¬ 
stream  temperatures  with  time.  It  is  this  record  that  provides  the  pri¬ 
mary  data  for  the  evaluation  of  heat  transfer  data.  The  precise  relation 
of  this  record  to  the  Colburn  j  factor  is  explained  in  Appendix  A.  The 
record  is  made  on  a  Leeds  and  Northrup  "Speedomax"  Temperature  Recorder. 

The  recorder  has  variable  chart  speeds,  which  permits  maintenance  of  accept¬ 
able  cooling  curves  with  regard  to  slope  determination.  The  Speedomax  also 
has  a  range  of  0-20  mv.  For  this  experiment,  the  recorder  is  precalibrated 
on  a  3  mv.  scale. 

There  is  one  addition  thermocouple  grid  that  is  not  used  in  the  experi¬ 
ment,  but  was  used  in  checking  out  the  apparatus.  This  grid  would  take  the 
place  of  T^  and  its  purpose  was  to  determine  the  temperature  distribution 
across  the  cross  section  in  both  the  vertical  and  horizontal  directions. 

It  was  placed  in  a  special  holder  that  also  held  a  velocity  probe  (See 
Figure  22)  and  was  lined  with  balsa  wood  to  provide  a  smooth  and  continu¬ 
ous  internal  flow  passage. 

PRESSURE  MEASURING  SYSTEM;  (See  Figure  1) 

Pressure  taps  are  located  upstream  and  downstream  of  the  orifice  and 
upstream  and  downstream  of  the  matrix.  Each  tap  is  connected  by  Saran 
tubing  to  an  appropriate  manometer  or  draft  gage.  The  following  instru¬ 
ments  were  used: 

1)  Type  "C"  micromanometer,  E.  Vernon  Hill  &  Co.  0-1.25" 

2)  Ellison  Differential  Direct  Draft  Gage,  0-6"^  0-8" 

3)  Ellison  Inclined  Draft  Gage,  0-6" 

4)  Meriam  Instrument  Co.  0-60". 
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Any  one  or  combination  of  these  instruments  can  be  used  to  measure 
the  differential  pressures  of  the  orifice  and  the  matrix.  Frequent  cross 
checks  were  made  to  ensure  reliable  operation. 

FLUID  HEATER  SYSTEM;  (See  Figure  1) 

The  air  heater  section  consists  of  28  nichrome  .0031"  diameter  wire 
heaters.  The  heater  system  is  designed  to  heat  1000  pounds  of  air  per 
hour,  20°F  above  ambient  temperature.  Nichrome  was  chosen  because  of 
its  high  resistivity,  low  thermal  conductivity  and  low  specific  heat. 

The  nichrome  wire  is  wound  on  a  bakelite  frame,  two  heaters  to  a 
frame.  The  wires  are  1/32"  apart  and  there  are  50  and  52  wires  respective¬ 
ly,  to  each  heater  on  a  frame.  The  heaters  on  each  frame  are  connected  in 
parallel  and  then  connected  through  a  switch,  to  a  variable  voltage  bus. 
There  are  two  frames  to  a  plastic  holder,  as  shown  in  Figures  23  and  24. 

In  these  figures,  the  heater  frame  and  its  electrical  connections  can  be 
clearly  seen. 

A  schematic  wiring  diagram  for  one  frame  is  shown  below.  All  14 
frames  are  wired  in  this  manner,  putting  all  28  heaters  in  parallel  with 
the  capability  of  switching  on  or  off  two  at  a  time. 
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The  28  heaters  arranged  in  this  way,  provide  a  variety  of  flexi¬ 
bility.  The  number  of  heaters  in  use,  should  decrease  as  the  flow  rate 
decreases.  Voltage  variation  can  be  obtained  by  a  GENEML  I?ADIO  COMPANY 
"VARIAC,  0-240'^,  8a,  50-60  cycles.  It  is  recommended  that  a  voltage  of 
220^^  or  in  this  vicinity,  be  maintained,  at  all  flow  rates,  and  the 
number  of  heaters  used,  decreased  until  the  number  of  heaters  is  four. 

At  this  point,  the  voltage  should  be  decreased,  as  necessary,  to  obtain 
a  20°  temperature  rise,  the  number  of  heaters  remaining  constant.  A 
note  of  caution;  Do  not  energize  a  heater,  unless  the  air  is  flowing, 
otherwise,  the  wires  will  sag  due  to  thermal  expansion  and  short  circuit¬ 
ing  will  occur. 

MATRIX  HOLDER  AND  TEST  SECTION; 

The  matrix  holder  and  test  section  casing  are  made  from  polyethylene 
plastic.  This  material  was  chosen  because  of  its  low  heat  capacity. 

The  test  section  casing  holds  thermocouples  T2  and  T^  permanently 
mounted  upstream  and  downstream  respectively.  It  also  contains  the  up¬ 
stream  and  downstream  static  pressure  taps  used  in  determining  the  fric¬ 
tion  factor.  The  parts  making  up  the  casing  were  machined  to  close  toler¬ 
ances,  to  ensure  a  tight  fitting  matrix  holder,  and  a  minimum  of  air  leaks. 

The  parts  comprising  the  matrix  holder,  as  implied  above,  were  machin¬ 
ed  to  close  tolerances  to  ensure  a  tight  fit.  Great  care  was  taken  to 
ensure  smooth  alighment  of  the  flow  channel  through  the  heaters  and  the 
test  section.  The  size  of  the  matrix  holder  was  arbitrarily  chosen. 

The  flow  channel  dictated  the  width,  and  the  length  is  such  that  matrices 
of  various  flow  lengths,  up  to  3",  may  be  tested.  The  flow  channel  is 
3-3/16”  X  3-3/16”. 

The  matrix  holder  has  one  special  provision  and  that  is  for  the  movable 
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thermocouple  A  small  hole  is  drilled  in  the  front,  and  the  wires 

are  led  through.  The  matrices  are  placed  in  the  holder  and  held  in 
place  by  bits  of  Styrofoam  plastic  cut  to  the  size  necessary  to  hold 
the  matrix  square  in  the  flow  channel. 

That  completes  a  description  of  the  various  systems  involved  in 
the  apparatus,  but  there  is  one  other  item  that  bears  mentioning.  Just 
upstream  of  the  heater  section,  is  a  section  provided  for  screens  to 
straighten  the  velocity  profile.  Provisions  were  made  for  1  to  4  screens. 
One  60  mesh  stainless  steel  screen  was  found  to  be  sufficient,  to  straight¬ 
en  the  flow  at  the  flow  rates  used  (40-1000  Ib/hr.) 
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APPENDIX  C 


CALCULATION  OF  GEOMETRIC  CONSTANTS 
FOR  PERFORATED  MATERIAL 

As  demonstrated  in  Appendix  A,  the  surface  geometry  has  a  strong 
influence  on  the  heat  transfer  and  friction  factor  curves.  Specifically, 
the  effects  of  p  and  ^  were  examined. 

In  using  perforated  materials,  the  presence  of  the  slots,  require 
modification  of  these  two  parameters.  It  is  the  intent  of  this  Appendix, 
to  demonstrate  the  calculation  procedure  utilized  to  determine  the  geo¬ 
metric  parameters,  on  which  the  graphical  and  tabular  results  for  the 
matrices  constructed  from  perforated  material  are  based. 

The  material  used  in  the  perforated  plate  matrices,  was  a  pure  nickel 
type  160/40TV,  slotted  opening  plate  material,  .0022  inches  thick.  This 
material  is  a  product  of  Perforated  Products,  Inc. 

The  basic  geometry  of  the  plate  is: 


^1— 


-ir/ff/K  "/////A  m/K 


li 

"T 


/////A 


[4.  .03BtM*U  .oi'Uj 


Using  the  average  between  the  top  and  bottom  of  the  tapered  slots, 
consider  an  average  slot  pattern  as: 
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CROSS  HATCHED  AREA  A  ,  =  (.063)  (.016)  =  .001008  cm^ 
cn 

AREA  OF  SLOT  IN  A^^  =  (. 040)  (. 0020)  =  .00008  cin^ 

Now  define  Y  =  Aslot  =  AREA  REDUCTION  DUE  TO  SLOTS 


.00008 

.001008 


7.937. 


Therefore,  plate  Porosity  (p^)  =  (1  -  Y  )  =  .9207 


Now  define 


S 


CROSS  SECTION  AREA  OF  SLOTS 

CROSS  SECTION  AREA 


CONDUCTION  AREA 
REDUCTION  DUE  TO 
SLOTS 


6- 


(.040) (.0022) 

(.063)(.0022) 


=  63.57o 


The  following  calculations  are  for  the  20°  skew  perforated  plate 


matrix: 


DIMENSIONS  OF  PLATE: 

4.285"  X  1.98"  X  .0022' 

NUMBER  OF  PLATES: 

142 

DIMENSIONS  OF  MATRIX: 

3.255"  X  3.248"  x  1.98' 

WEIGHT  OF  MATRIX: 

275.5735  grams 

MATERIAL  CONSTANTS; 

^  =  555  lb/^^3 

ks=  38.7  BTU/^^„^ 

Cs=  -10^5  BTU/^j^op 

PLANE  SURFACE  AREA  =  (1 . 98)  (2)  ^  16.733  ft^ 

144 

TOTAL  HEAT  TRANSFER  AREA  (A)  =  PLANE  SURFACE  AREA  -  AREA  OF  SLOTS 
=  (1-  ^(*  )  SURFACE  AREA 
=  (.9207)(16.733) 


FRONTAL  AREA  (A^^) 


A  =  15.41  ft 


3.255  X  3.248  =  10.572 


,  2 
in 
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FREE  FLOW  AREA  =  FRONTAL  AREA  -  CROSS  SECTIONAL  AREA  OF  METAL 

2  2 
=  10.572  in  -  4.285  x  .0022  x  142  =  9.233  in 


A  =  .06416  ft^ 
c 


CONDUCTION  AREA  (A^)  =  (142)  (4. 285) (.0022)  (1- ^ ) 

=  (142)(4.285)(.0022)(.375)  =  .50198  in^ 


A  =  .003486  ft^ 

s  _ 


MAXTIX  VOLUME  (V^)  =  3.255"  x  3.248"  x  1.98"  =  20.933  in^ 

MATERIAL  VOLUME  (V^)  =  142x4. 285"xl . 98"x. 0022"  =  2.650  in^ 

L'EIGHT  OF  MATERIAL  (W^)  =  275.5735  gms./453.6  gm/lb  =  .6075  lb. 

MATERIAL  VOLUME  (T)  =  W  /  «  =  .6075  lb/555  Ib/ft^  =  .001095 

s  s  Ps 

The  porosity  of  a  matrix  has  been  defined  as: 


MATRIX  FLOW  VOID  VOLUME 

MATRIX  VOLUME  "  m  ~  s 

Vm 


Vs 

V 

m 


P 


.001534 

.01211 


.8733 


P 


.  .001095  _ 

V  ^  .01211  "  ^ 


.09042  =  .90958 


The  compactness  of  a  matrix  has  been  defined  as 


/S 


TOTAL  HEAT  TRANS FE R  AREA 
MATRIX  VOLUME 


Here, 


we  will  redefine 


to  include  the  effect  of  area  reduction: 


(PLANE  SURFACE  AREA)(1-  V  ) 

MATRIX  VOLUME 
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For  '^  =  0 


(5  = 


(16-733)(1) 

.01211 


For 


^  =  .0793 


(16.733)C.9207) 

.01211 


1272.18 


The  hydraulic  diameter  is  conventionally  defined  as; 


4  FLOW  CROSS  SECTION  AREA 
^  WETTED  PERIMETER 

Appendix  A  shows  that  hydraulic  diameter  can  also  be  expressed  as: 


Dh  =  4 

To  be  consistent  with  the  conventional  definition,  for  the  case  of 
a  perforated  material,  ^  is  used  vice  ^  ,  since  the  slots  have  no 

effect  on  the  hydraulic  diameter. 

Finally  then. 


.8733 

1381.75 


.002528  ft 
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APPENDIX  D 


SAMPLE  CALCULATIONS 


The  following  sample  calculations  are  based  on  recorded  data  from 


Run  43  for  20“  skew  matrix  constructed  from  perforated  nickel  sheets. 
This  is  chosen  primarily  because  the  basic  geometric  parameters  are  com¬ 


puted  in  Appendix  C. 


RECORDED  DATA: 


D 

2.310” 

o 

> 

o 

= 

2.61  in. 

H2O 

AP 

3.96  in. 

H.,0 

m 

2 

P 

_ 

4.69  in. 

H_0 

s 

2 

T 

= 

68.0  “F 

MATERIAL  CONSTANTS: 

c^  =  .1065  BTU/lb  “F 

k^  =  38.7  BTU/hr  ft“F 
PRECALCULATED  GEOMETRIC  PARAMETERS; 

A  =  15.41  ft^ 

A^  =  .0734  ft^ 
fr 

A  =  .06416  ft^ 
c 

r  =  .003486  ft^ 
s 


P  ^  =  29.963  in.  Hg 

atm  ® 

CHART  SPEED  =4.0  in/ sec 
W  r.6075  lb 


(Re:  Appendix  C) 
lO  =  .8733 

|3  =  1272.18  ft^/ft^ 
r^  =  .000632  ft 


L  =  1.98  inches  (  —  denotes  corrected  for  slot  effects) 
DETERMINATION  OF  MASS  FLOW:  as  outlined  in  Reference  (7) 

m  =  359  K  D  ^  F  Y  J  A  P  T” 
o  a  ^  o 

where;  Y  *  1.0  Fig.  40b  Ref.  (7) 

Fg  «  1.0  Fig.  38  Ref.  (7) 
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in 


(359)(2.31)^  K  (2.61)(  (.07522) 


m  =  849.40  K 


assume  K  =  .749 


Then  m  =  (849.40) (.749)  =  636.21  Ib/hr 

Reynold's  number  for  the  pipe  N  = 

Rp 


=  7.35  X  10 


From  Table  5,  Ref.  (7)  0.75,  =  7.35  x  10^,  K  =  .749 


Therefore 


m  =  636.21  ^^/hr 


MATRIX  REYNOLD'S  NUMBERS 

There  are  actually  two  Reynold's  numbers  to  be  calculated  for  the 
matrix.  One  is  for  isothermal  flow  and  is  associated  with  the  .friction 
factor,  the  other  is  associated  with  heat  transfer.  They  differ  only  in 
the  temperature  at  which  the  absolute  viscosity  ^  is  evaluated.  ^ 

(friction)  is  evaluated  at  ambient  temperature  (T^) .  ^  (heat  transfer) 

is  evaluated  at  T^  +  10®F  (an  estimate  of  the  average  bulk  fluid  tempera¬ 
ture)  . 


(^)(^^)  ft (636. 21^^,  hr 


=  618.58 


(15. 41ft^) (.04405  Ib/ft  hr) 


CONDUCTION  PARAMETER 


k  A 


.00536 
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MAXIMUM  SLOPE 


From  Figure  26  the  slope  of  the  cooling  curve  was  determined  as 
described  in  Appendix  A.  The  slopes  of  three  curves  were  determined  and 
the  average  of  the  three  was  used  in  the  calculations. 


MAXIMUM  SLOPE  = 


^s  /  dT  \  _  ^s 


SLOPE  xCHART  SPEED 


^  d0/max 


(.6075  lb)(.1065  “^^/lb°F)  =  .06470'^^^/°F 


C  =  m  Cp  =  (636.21^^/hr)(.24®‘^^/lb°F) 


.65745 


SLOPE  =  ,16014  in  (See  Figure  26) 

Finally: 

MAXIMUM  SLOPE  =  (*10775)  4) 

COLBURN  j  FACTOR 

With  X  and  MAXIMUM  SLOPE  enter  Table  VIII  and  find  N  =  4.40. 

tu 

Linear  extrapolation  will  give  desired  accuracy. 


Now, 


J  =  r  '•pR 


2/3 


from  Appendix  A,  Eq.  (A-22) 


therefore,  j  =  *  15^41  (.796)  (4.40)  = 

FANNING  FRICTION  FACTOR 


.01458 


As  derived  in  Appendix  A  Equation  (A-16)  the  friction  factor  is: 

AP  =  3.96  in  H„0  x  5.204  ■  =  20.61  psf  =  .1431  psi 

m  2  in  H,0 
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AF  =  4.69  in  H„0  x  5.204  =  24.41  psf  =  .1695  psi 

s  4  In  H2O 

P  =  29.963  in  H  X  .4912  =  14.7178  psi 

atm  g  in 

P  =  P  -  P  =  14.7178  -  .1695  =  14.5483  psi 
1  atm  s 


14.4768  psi 


^  ^  (U.477)(144)  ^  lb  3 

rmean  RT  (53.3)(528) 

I  m 

m  636.21  Ib/hr  «  ,,  /,  2  ,^,lb,  ^.2 

^  ~  =  .064i6"'ft2~  =  =  2-754  /sec  ft 


.000632 
Vl  “  1.98/12 


=  .00383 


K  =  -.738 


K  =  .949  from  Chapter  IV  Reference  (8) 


(1  +  p  )  =  1  +  (.8733)  =  1.7626 

AP 


28. 


mean 


_  _  (64.4)(.07407)(20.61)=  12.9579 

,2  ■  7.587 


AP  r,  2. 

_ni _  (1+P  ) 


■  a-7«)  =  -oi^w 

14.4768  psi  ' 


f  =  (12.9579  -  .01740  -  . 211) (.00383)  =  (12.7295) (.00383) 

=  .04859 

All  thermal  properties  (  ,  Npj^,kg,  and  c^)  jere  obtained  from 

References  (12)  and  (13).  The  temperature  at  which  the  properties  were 
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evaluated  was  the  ambient  temperature  with  the  exception  of  the  absolute 
viscosity  for  heat  transfer  (  )• 
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APPENDIX  E 


A  FORTRAN  COMPUTER  PROGRAM 
FOR  DATA  REDUCTION 

A  short  computer  program  was  devised.  Figure  27,  so  that  from  basic 
experimental  data,  the  end  result  is: 

1)  Conduction  parameter  (A) 

2)  Maximum  slope 

3)  Reynold's  number  for  heat  transfer 

4)  Reynold's  number  for  friction 

5)  Fanning  friction  factor 

It  was  anticipated  that  this  program  along  with  a  subroutine  would 
ultimately  produce  N^^^  and  subsequently  the  Colburn  j  factor. 

Time  did  not  permit  the  adaption  of  Moreland's  work  (11)  to  a  sub¬ 
routine  but  the  program  is  still  useful.  It's  value  when  data  for  a 
single  matrix  is  to  be  reduced  is  nebulous.  However,  the  reduction  of 
data  for  multiple  matrices  can  be  accomplished  with  a  considerable  savings 
of  time. 

The  following  is  a  list  of  inputs  necessary  in  Fortran  and  English 
language  along  with  the  required  dimensions: 

CARD  NO. 


XML 

Matrix  length 

ft 

AS 

Solid  conduction  area 

ft^ 

RH 

Hydraulic  radius 

ft 

POR 

Porosity 

- 

AFL 

Flow  cross  section  area 

ft^ 

AHT 

Matrix  heat  transfer  area 

ft^ 

WM 

Matrix  weight 

lbs 
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CARD  NC. 


40  EK  Exit  flow  coefficient 

CK  Entrance  flow  coefficient 

SK  Solid  thermal  conductivity 

CM  Matrix  material  specific  heat 

N  Number  of  runs 

10  DO  Diameter  of  orifice  plate  opening 

DELPO  Pressure  drop  across  orifice 

GAM  Specific  weight  of  fluid 

ATMP  Atmospheric  pressure 

COR  Orifice  correction  K 

UH  Fluid  viscosity  for  heat  transfer 

UFR  Fluid  viscosity  for  friction 

10+2  CS  Chart  speed 

DELH  Pressure  drop  across  matrix 

HS  Static  pressure  matrix  inlet 

SLO  Slope  of  generalized  cooling  curve 

TEMPO  Orifice  temperature 

The  output  of  this  program  in  Fortran  and  English 
their  dimensions  is: 

EMDOT  Mass  flow  rate  (m) 

RNUMP  Reynold’s  number  for  the  pipe 

RNUMH  Reynold's  number  for  heat  transfer 

CAPFL  Thermal  capacitance  of  the  fluid  (C) 

CAP^  Thermal  capacitance  of  the  solid  (C^) 

CONDPAR  Conduction  parameter  (  X  ) 


BTU/hr“F  ft 
BTU/lb“F 


inches 
inches  H2O 
Ibs/ft^ 
inches  Hg 

Ib/hr  ft 
Ib/hr  ft 
sec/ inch 
inches  1120 
inches  H2O 
(Inches)  ^ 

“F 

language  along  with 

Ib/hr 


BTU/sec  “F 
BTU/  “F 
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SLOMAX 


Maximum  slope  of  generalized  cooling  curve 


G 


Mass  flow  velocity 


Ib/hr  ft' 


2 


RNUMFR 


Reynold's  number  for  friction 


FFR 


Fanning  Friction  factor 


Reduction  of  data  is  not  complete  since  j  remains  to  be  found. 
With  SLOMAX  and  CONDPAR,  enter  Table  VIII,  and  find  N  .  The  Colburn  j 
factor  follows  from  equation  (A-22),  Appendix  A. 

The  program  is  flexible  in  that  an  unlimited  number  of  matrices 
may  be  reduced  by  changing  the  first  card  of  the  program  appropriately. 
(L  =  the  number  of  matrices  to  be  reduced.) 
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